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ABSTRACT 


AMPLIFICATION  OF  SOUND  BY  GAS-PHASE  REACTIONS 

DETSCH,  RICHARD  MARK.  B.S.,  Pennsylvania  State  University,  1977.  M.S., 

Western  Kentucky  University,  1979.  Ph.D.,  University  of  Mississippi, 
1983.  Dissertation  directed  by  Professor  Henry  E.  Bass. 

A  four  year  study  of  sound  propagation  in  chemically  reacting  mixtures 
has  led  to  experimental  observation  of  sound  amplification.  Photo- 
initiated  Cl2-H2-inert  gas  reactions  provided  the  energy  for  the  amplifica¬ 
tion  observed.  Amplification  experiments  were  conducted  in  two  modes  which 
we  termed  pulsed  and  cw.  For  the  pulsed  experiments,  high  intensity  UV- 
flash  lamps  dissociated  molecular  chlorine  at  the  start  of  the  experiments 
and  the  Cl2-H2-Ar  reactions  quickly  went  to  completion;  amplification  was 
observed  for  about  6  msec.  The  measured  gain  was  1.8  at  2  kHz  compared  to 
a  gain  of  4  predicted  by  theory.  The  source  of  this  difference  has  not 
been  resolved  completely. 

In  cw  experiments,  low  intensity  UV-fluarescent  lamps  dissociated 
molecular  chlorine  continuously  during  the  experiments,  this  extended  the 
time  over  which  amplification  was  observed  to  about  800  msec.  Experiments 
were  conducted  for  three  C12-H2-SF6  mixtures  and  four  frequencies,  1.0, 

2.5,  4.0,  and  6.5  kHz.  Computer  calculations  correctly  predicted  the 
frequency  dependence  of  amplification  and  gave  amplification  values  which 
were  in  fair  agreement  with  observation. 


Attenuation  measurements  were  conducted  for  several  CI2-SF6  and  H2- 
SFe  mixtures.  These  data  were  used  to  calculate  vibrational  relaxation 
times  for  SFs  in  collisions  with  Cl2,  H2,  and  SF6  molecules.  The  vibra¬ 
tional  relaxation  times  were  used  to  predict  vibrational  relaxation  atten¬ 
uation  in  CI2-H2-SF6  mixtures  used  in  cw  amplification  experiments. 
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CHAPTER  1 


INTRODUCTION 

The  interaction  between  nonequilibrium  chemical  reactions  and  sound 
waves  propagating  in  the  reacting  mixture  affect  a  variety  of  observable 
phenomena,  e.g.  combustion  instability  in  jet  and  rocket  engines,  struc¬ 
ture  of  detonation  waves,  and  turbulence  in  reacting  flows.  Studies  which 
isolate  the  effects  of  the  interaction  are  quite  limited.  There  is  reason 
to  believe1  that  this  interaction  might  lead  to  several  new  phenomena  or 
explanations  of  unexplained  observations.  Among  these  are  acoustic  stimu¬ 
lation  of  chemical  oscillations  and  chemical  instability,  amplification 
of  sound,  changes  in  sound  speed  and  frequency  during  the  course  of  the 
reaction,  and  sound  induced  changes  in  reaction  rates.  In  this  disserta¬ 
tion,  amplification  of  sound  will  be  of  primary  interest.  While  some  phe¬ 
nomena  may  be  observable  at  near  equilibrium  conditions2-4 ,  extreme  non- 
equilibrium  conditions  are  required  for  measurable  sound  amplif ication1 . 
There  have  been  a  number  of  theoretical  studies  of  sound  amplification  in 
irreversibly  reacting  mediums5-13.  Experimental  verification  of  this 
phenomenon,  however,  is  limited6.  Related  mechanisms  for  producing  sound 
amplification  in  a  gas  include  ionization  of  a  gas14,  pressure  dependence 
of  infrared  absorption15,  and  selective  excitation  of  internal  modes16. 
Experimental  verification  of  these  mechanisms  is  also  limited. 

Devices  which  amplify  sound  directly  will  be  refered  to  in  the  fol¬ 
lowing  as  SACERS  (Sound  Amplification  by  Controlled  Excitation  Reactions) . 


The  purpose  of  this  investigation  was  to  observe  sound  amplification  in  a 
chemically  reacting  gas  mixture  (a  chemical  SACER) .  The  H2-Cl2-buf fer 
gas  photoinduced  chemical  reaction  was  chosen  for  several  reasons:  the 
gases  could  be  mixed  without  reacting  and  then  UV-irradiated  to  initiate 
the  reactions  at  a  controlled  time;  most  of  the  reaction  rates,  activa¬ 
tion  energies,  etc.  were  known;  and  several  theoretical  studies  had  been 

done  on  this  reaction  scheme5-8.  Acoustic  frequencies  in  the  range  of 

1 

1-lOkHz  were  chosen  because  lower  frequencies,  having  periods  longer  than 
the  relaxation  times  of  the  reactions,  introduce  theoretical  complications 
absorption  calculations  in  this  range  have  been  done  previously  in  this 
lab;  and  higher  frequencies  have  increased  attenuation  making  it  more 
difficult  to  observe  amplification. 

Two  SACERs  were  constructed  for  this  investigation.  The  first  con¬ 
sisted  of  a  Cl2-H2-Ar  gas  filled  tube  surrounded  by  high  intensity  flash 
lamps.  The  tube  was  transparent  to  UV-radiation  except  for  small  unillu¬ 
minated  regions  at  each  end.  The  flash  lamps  completely  dissociated  Cl2 
in  the  illuminated  region  causing  a  pressure  increase  which  expanded  into 
the  unilluminated  regions.  This  expansion  led  to  pulses  which  propagated 
down  the  tube.  A  microphone  recorded  the  amplitude  of  these  pulses  at 
various  propagation  distances.  Electronic  filters  were  used  to  study 
individual  frequencies  found  in  the  pulses.  This  SACER  is  described  as 
operating  in  a  pulsed  mode  because,  after  flashing,  the  reactions  lasted 
only  about  10  msec  (the  time  required  to  use  all  the  H2). 

The  second  SACER  used  low  intensity  flourescent  lamps  to  dissociate 
Cl2.  Since  only  a  small  amount  of  Cl2  was  dissociated  at  any  given  time. 


the  reactions  continued  much  longer  (approx.  1  sec);  we  refer  to  this  as 
the  cw-SACER  (continuous  wave) .  Another  reason  for  the  slower  reaction 
rate  was  the  substitution  of  SFg  for  Ar.  SF6  has  a  larger  specific  heat 
which  slowed  the  temperature  increase  of  the  gas.  By  keeping  the  transla¬ 
tional  temperature  low  the  reactions  evolved  less  rapidly.  In  this  system 
there  were  no  unilluminated  regions;  tone  bursts  were  generated  by  a  send¬ 
ing  transducer.  A  receiving  transducer  measured  the  amplitude  of  the  tone 
burst  after  propagating  various  distances  through  the  reacting  gas. 

A  necessary  part  of  this  study  was  measurement  of  sound  attenuation 
in  various  SF6-CI2  and  SF6-H2  mixtures.  From  these  measurements  vibra¬ 
tional  relaxation  times  of  SFg  could  be  predicted  for  the  CI2-H2-SF6  mix¬ 
tures  which  made  it  possible  to  calculate  the  attenuation  due  to  vibra¬ 
tional  relaxation  in  cw-SACER  experiments. 

Chapter  II  presents  a  mathematical  and  physical  description  of  sound 
amplification  in  chemically  reacting  systems.  Also  described  is  the 
theory  of  sound  attenuation,  including  that  due  to  vibrational  relaxation. 
The  theory  is  not  original  to  this  work  and  is  presented  only  for  the 
benefit  of  the  reader  unfamiliar  with  previous  work  in  the  field.  This 
study  did  require  that  the  gain  calculation  of  Gilbert5*6  be  meshed  with 
the  attenuation  calculations  based  on  work  by  Shields17.  That  step  will 
be  described  in  more  detail.  The  primary  thrust  of  this  study  was  collec¬ 
tion  and  analysis  of  experimental  data.  Most  effort  was  expended  in  con¬ 
struction  of  the  experimental  apparatus  described  in  Chapter  III.  Chapter 

III  includes  a  description  of  the  acoustic  waveguides,  the  UV-radiation 
sources,  and  the  computer  system  which  controlled  the  cw-SACER.  Chapter 

IV  gives  experimental  results.  Conclusions  are  discussed  in  Chapter  V. 


CHAPTER  II 


THEORY 

Introduction 

Two  conditions  can  generate  sound  amplification  in  gas-phase,  chemi¬ 
cally  reacting  systems.  The  first  occurs  when  the  reaction  has  an  imbal¬ 
ance  of  stoichiometric  coefficients  and  a  pressure  dependent  reaction  rate. 
Sound  modulates  the  pressure  causing  spacial  fluctuations  in  the  reaction 
rate.  The  rate  fluctuations  are  accompanied  by  fluctuations  in  total  atom 
concentration  and  hence  secondary  pressure  fluctuations.  As  an  example 
consider  the  reaction 

A  — - — »-  B  +  C  (1) 

where  the  reaction  rate,  R(p)  increases  monotonically  with  pressure.  As 
the  reaction  progresses,  the  overall  pressure  is  increasing,  however,  the 
rate  of  pressure  increase  is  modulated  by  the  acoustic  signal.  The  reac¬ 
tion  rate  is  larger  in  compressed  regions  causing  a  relative  increase  in 
pressure  due  to  the  imbalance  of  stoichiometric  coefficients.  The  increase 
in  pressure  caused  by  the  sound  results  in  a  feedback  mechanism  yielding 
a  further  increase  in  pressure.  A  similar  feedback  mechanism  occurs  in 
rarefied  regions  thus  the  acoustic  signal  is  amplified.  Reversing  the 
pressure  dependence  of  the  reaction  rate  or  the  stoichiometric  imbalance 
results  in  negative  feedback. 


The  second  amplification  mechanism  occurs  when  the  reaction  is  exo¬ 
thermic  and  has  a  reaction  rate  that  increases  with  temperature.  Under 
adiabatic  conditions,  sound  modulates  the  temperature  and  thus  the  reac¬ 
tion  rate.  In  regions  of  positive  temperature  change,  reaction  rate  and 
heat  generation  are  increased  causing  a  secondary  temperature  increase. 
Reversing  the  temperature  dependence  of  the  reaction  or  having  an  endo¬ 
thermic  reaction  causes  negative  feedback. 

These  two  feedback  mechanisms  can  be  present  simultaneously.  They 
can  also  be  coupled,  i.e.  a  pressure  dependent  reaction  which  is  exother¬ 
mic.  Note  that  normal  attenuation  mechanisms  are  still  present  and  can 
cause  an  absolute  decrease  in  sound  intensity,  with  propagation  distance, 
even  when  positive  feedback  is  present. 

In  the  following  the  general  formalism  of  non-equilibrium  chemical 
dynamics  necessary  to  treat  sound  propagation  in  a  chemically  reacting 
mixture  is  described.  This  treatment  is  available  in  the  cited  literature 
so  only  the  basic  physics,  major  conclusions,  and  details  useful  to  future 
students  will  be  included.  In  practice,  all  computations  were  performed 
using  a  program  written  by  Professor  Gilbert13.  Next,  mechanisms  for 
attenuation  are  identified  and  described  only  to  the  detail  required  for 
consistency.  A  program  written  by  Professor  Shields17  (see  appendix  A) 
was  used  for  this  calculation  with  a  modification  to  allow  viscosity  to  be 
computed  as  a  function  of  temperature  for  a  ternary  mixture. 


General  Formalism  for  Non-Equilibrium  Thermodynamics 


Equations  needed  to  describe  chemical  amplification  are  presented  in 
this  section.  The  following  is  taken,  in  large  part,  from  Reference  13. 
Consider  the  reaction 


2H2  +  02  -*  2H20 


which  is  the  form 


2A1  +  A2  -+•  2Aj . 


If  this  is  the  k1-*1  reaction  being  considered,  the  actual  stoichiometric 
coefficients,  nuk*  are  ^ik  "  "2,  n2fc  ■  -1,  and  ri3k  *  +2.  Eq.  (3)  can  be 
written  in  summation  form  as 


An«.kA2,  ”  0 

l 

8 

where  £  indicates  a  summation  over  all  species.  In  terms  of  grams, 
instead  of  moles,  eq.  (3)  is 


4Aj  +  32A2  -*■  36A3 


where  Aj  represents  one  gram  of  species  j.  Reducing  eq.  (5)  gives 
A3  -  Ai  -  ■—  A2  "0 


IvikAi  “  0 


'Vnrr? 


7 

where  vik  is  the  mass  weighted  stoichiometric  coefficient  for  the  ifch 
species  and  reaction.  Let  the  progress  variable  of  the  kth  reaction, 
Xk»  be  defined  as  the  number  of  grams  of  reactant  which  have  reacted  per 
gram  of  original  reactant.  Then 


dfcmi  -  mvikdXk 


(8) 


where  m^  is  the  total  mass  of  species  "i",  m  the  total  mass,  and  d^m^  the 
change  in  m^  caused  by  the  kc^  reaction.  The  total  change  in  mi,  dmi, 
caused  by  all  the  reactions  is 


r  r 

dmi  "  I^kmi  -  mJvikdXk* 
k  k 


(9) 


Dividing  eq.  (9)  by  the  total  volume  gives 


dPi  -  P^VikdXk 
k 


(10] 


where  Pi  is  the  mass  density  species  "i",  and  p  the  mass  density. 
If  one  considers  only  forward  reactions  then,  for  our  example 


pdXj/dt  -  p£j  -  kjpfp2 


(11] 


where  is  the  degree  of  advancement  of  the  kfch  reaction,  =dXj/dt,  and 
kj  is  the  rate  coefficient  of  the  jt*1  reaction  in  the  forward  direction. 
In  general  form, 


-  p”  kjpi  '1p2 


(12] 


or,  in  terms  of  mass  fraction,  c j , 


€j  -  p-ikjCrnic2_T,2p"ni‘n2-p(8j“1>kjcrnic2_n2 


(13 


where  gj  is  the  molecularity  of  the  reaction.  For  reactions  of  any 


order 


-  p(ej-»kJJ(c1)-'’i  <u) 

1  s' 

where  is  the  mass  fraction  of  the  ifch  reactant,  and  £  the  sum  over  all 

reactants . 

Consider  any  extensive  variable,  6,  i.e.  energy  or  mass.  Let  G  be 
specific  6,  6  per  unit  mass.  Then 


0  -  J pGdv 


where  the  integral  extends  over  the  entire  system;  and 


0  ■  d/dt  JpGdv. 


If  the  volume  remains  constant,  the  time  derivative  may  be  taken  inside 
the  integral  giving 

0  «  VJ~  (pG)dv.  (17) 

An  alternate  expression  for  6  can  be  derived.  If  a  flow  of  0  is  present 
at  the  surface  of  a  volume,  then  the  resulting  change  9  with  time  is 
-aJjG*d8,  where  Jq  is  the  current  density  for  the  extensive  variable  9, 
pG9;  and  v  is  the  center-of-mass  velocity.  If  there  are  internal  sources 
of  0,  denoted  $q,  the  resulting  change  in  9  with  time  will  be  vJ$gdv. 

Thus, 


9  "  -apG*dS  +  vJ<Mv. 


Equating  the  two  expressions  for  9  gives 


where  Jgx>m,x  is  the  x-component  of  momentum  current  density,  »  pvxv;  Xx 
the  x-component  of  external  force  per  unit  volume;  and  a  the  stress  tensor 
Combining  eq.  (26)  with  its  y-  and  z-  counterparts  gives 

-^•(pv)  +  V •  (pvv)  -  pX  +  V*a  (27) 

and,  using  conservation  of  mass, 

pv  *  pX  +  V •£.  (28) 

Assuming  no  external  forces  and  using  the  stress  tensor  for  homogeneous, 
isotropic  fluids,  eq.  (28)  can  be  written 

pv  ■  -Vp  +  nV2v  +  (jn  +  4>)V(V*v)  (29) 

where  n  is  the  coefficient  of  shear  viscosity,  and  $  the  coefficient  of 
bulk  viscosity.  Applying  eq.  (20)  to  entropy  gives 

pS  -  $/ T  -V*JS  (30) 

where  S  is  the  specific  entropy;  $/T  the  source  potential  of  entropy;  and 
js  the  entropy  flux  caused  by  diffusion  and  heat  flow.  4>  is  caused  by 
viscous  flow,  chemical  reactions,  isothermal  diffusion  and  heat  flow. 
Fitts18  calculated  $,  and  by  using  conversation  of  mass  and  eq.  (30), 
obtained  the  following 

r 

pK"Cvp  ■  pCp7*v  -  £pSj(gAHj  -  pCpAVj)  (31) 

j 

and 

r 

pK'Cy  T  -  -8T  V*v  -  ^p?j (K'AHj-STAVj ) 


(32) 


“»«'  4%  -  [  v^Si.  iVj  -  J  vy^,  K'  -  ^)T,H,  and  6  -  ±(|f>p,m. 
Eqs.  (31)  and  (32)  were  derived  using  the  following  assumptions: 


1.  adiabatic  conditions ,  7*q  *  0 

2.  no  diffusiona  ji  ■  0 

3.  homogeneous,  isotropic  fluid  with  bulk  and  shear 
viscosities  independent  of  gradients. 


The  necessary  acoustic  relationships  follow  from  the  conservation 
equations.  It  will  be  assumed  that  no  chemical  reactions  are  occurring 
and  that,  except  for  thermal  conduction,  no  attenuation  mechanisms  are 
present.  Attenuation  will  be  included  in  later  sections.  With  these 
assumptions,  eq.  (31)  becomes 


p  *  ~(cp/cvK')V*v 


For  plane  waves  propagating  in  the  x-direction  in  an  ideal  gas,  eq.  (33) 
reduces  to 


P  -  -(cp/cvK')p3vx/3x 


where  the  difference  between  the  substantial  derivative,  d/dt,  and  the 
time  derivative,  3/3t,  is  second  order  and  has  been  ignored.  With  the 
same  assumptions  as  those  used  to  obtain  eq.  (34),  eq.  (32)  becomes 


■jTr  T  “  -p(pcv)  x3vx/3x. 


eq.  (23)  becomes 


and  eq.  (28)  becomes 
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p7t  ^  “  (V‘2>X‘ 

Writing  the  stress  tensor  as 


(37) 


£ij 


■  -P^ij* 


eq.  (37)  simplifies  to 


.1.  8 
“'p'  3x  p* 


(38) 


(39) 


So  far,  we  have  developed  acoustic  and  chemical  relationships  between 
pressure,  center-of-mass  velocity,  temperature,  and  mass  fractions.  These 
expressions  have  been  combined  by  Gilbert  et  al.5  as  follows. 

Let  \|;(x,t)  be  the  space-  and  time -dependent  vector  of  variables 
describing  the  system.  The  components  of  this  vector  are  (p,  v,  T,  cj.) . 
Equations  describing  the  system  [eqs.  (25),  (29),  (31),  (32),  (34),  (35), 
(36) ,  (39)  and  (40) ]  can  be  combined  into  one  tensor  equation 


*  -  FMO  +  BfipJ  JL  f  (40) 

where  F[iJ>]  is  a  nonlinear  vector  function  describing  the  kinetics,  and 
B[i|»]  a  matrix  function  of  ^  describing  the  acoustics. 

The  quantity,  ip(x,t)  can  be  written  as  the  sum  of  ip°(t),  the  part 
resulting  from  homogeneous  evolution,  and  5\{/(x,t),  a  small  quantity  result 
ing  from  acoustic  perturbations.  Homogeneous  evolution  is  described  by 


i^°(t)  -  F[^»°] 


(41) 


or,  without  tensor  notation,  by 


P°  -  -p8^)-1  ic°t(T#)-1AHj  -p°CpAV|] 

Z 


v°  -  0 


T°  ■  -p8(cv)“l|c°[ (p°)-1AHj  -  AV£] 
Z 


Ivi zZ° • 


Inhomogeneous  evolution  is  described  by  eqs.  (25)  (29),  (31),  and  (32) 
These  four  equations  can  be  rewritten  as 


&  Ci  "  ci  + 

k 

3  _  -i  3  _  3  - 

at  vx  -  -P  3JP-Vx^vx 


3  d 

P  -  _vx  __  P  +  Cp  P  Cy-1-^-  Vx  -  pcy-^^kGk 


-i  3  _ 


■k T  ■  -**  k T  -  k  -  cv‘1|5jJj 

where  the  following  assumptions  and  definitions  have  been  made: 


rate  coefficients  have  an  Arrhenius  temperature 
dependence,  ki  -  f±  exp  (-E^/RT) ,  where  f±  is 
the  frequency  factor,  and  the  activation 
energy  for  reaction  "i"; 

the  molar,  constant-pressure  specific  heats  have 
the  form  c£  m  a±  +  gjT  +  nil”2,  where  a^,  6^, 
and  Hi  are  the  specific  heat  coefficients  for 
the  ith  species; 

Cp  ■  IcpCiM^-1 ,  where  Mg  is  the  molecular  weight 

of  species  "i". 

q±  5  °iT  +  JsSil'-mT-1 

Eg  S  E^/RT2-!”1 (gg-1) 
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J&  =  4%  ~  PAVj, 

Gjj,  =  T  —  pCpAVj^ 

=  molecular  weight  on  right-hand  side  of  reaction 
"4" 

s 

AH^CT)  -  CM^-HHfcCr  -  298K)  -  Ynit[qi(T  -  298K)-qi(T)  ] } 
Ai£  "  ||j“^4[-(5iir1iJl/ci)  -  (gt-l)Mi-1pRTp“1] 
ni£,  “  vi£M^/M£* 


Next,  let  the  independent  variables  be  composed  of  a  homogeneous  part 
plus  a  small  inhomogeneous  part,  i.e.  p  *  p°  +  6p.  Then  the  dependent 
variables  can  be  written  as  a  homogeneous  part  plus  a  Taylor  series  of 

Q 

inhomogeneous  variables,  i.e.  p  *  p°  +  (•|^-)“6p  +  (l=r)°(5T  +  £  (~-)°6ci 

n  i  dci 


+  (-r-*— )°<5vx,  where  higher  order  terms  have  been  neglected.  Using  this 

dVx 

notation  and  neglecting  higher  order  terms,  eqs.  (46)  -  (49)  become 

■&  5ci  "  5T  bikt&l  +  <«P0)_1I(g4-l)Uvi£  +  ^CjIvikAjk  (51 

k  l  j  k 

H  VX  -  0  (5 

^  6p  -  6p[-(c^)‘1ig)iqG|]  +  6T(-Cy)_1{p®f [G|5| 

+  (T°)-2(AH*-T  ^  AH*)°  -  P°AVJ(^)°}} 

s  r 

-  c$p0(9cv/3T)°6T  +  (-Cv)-l^6ck{p°^[G|  Ak* 

k  si 

+  CjdV|Mk~lp0(Cp/Y°-c^)  ]  }  -  c^-pcic^r1  (5! 

T  -  -(c^)_1{8T[ (XCj  C-^-  AH£)0-p8(T°)-1AVj  + J^Ep-(c^)2T°] 

+  fip[IJj(gA-l)(p0)_15j  +  ^5cit  ^Aj^Jj-(cy) 2CpMi-1T° ] }  (5: 

l  i  l 

where  the  homogeneous  parts  have  been  eliminated  using  eqs.  (41)  -  (45). 
Eqs.  (50)  -  (53)  can  be  written  in  matrix  notation  as 


■gt 


where  Q(^°)  is  a  matrix  whose  elements  are  found  by  solving  the  equation 
for  homogeneous  evolution,  (41).  A  numerical  solution6  of  eqs.  (41)  and 
(54),  written  by  Professor  Gilbert,  was  used  to  predict  sound  amplifica¬ 
tion  for  our  experiments.  The  solution  uses  the  Gear  method19  of  solving 
stiffly  coupled  differential  equations.  Another  study  of  sound  amplifica¬ 
tion  in  CI2-H2  gas-phase  chemically  reacting  mixtures  has  been  published 
by  Toong  et  al.7-9;  it  will  be  considered  briefly.  While  similar  to 
Gilbert's  solution,  Toong 's  derivation  neglects  viscosity  and  thermal  con¬ 
ductivity;  however,  this  simplification  allows  for  analytical  solutions. 
Thus,  clear  insights  concerning  acoustic-chemical  interactions  are  possi- 

a 

ble.  Toong  presents  a  useful  equation  for  comparing  sound  attenuation 
with  and  without  chemical  amplification 

Act  -  (2yT)”1dT/dt[m"+(EA/RT)  (y-1)  (Ti/T)+rl  (55) 

where  Tj,  is  the  initial  temperature  and  m'  the  order  of  the  reaction. 

Both  Gilbert  and  Toong 's  solutions  are  compared  to  our  experimental 
results . 

Attenuation 

In  order  to  predict  the  absolute  sound  amplitude  in  a  reacting  mixture 
the  standard  attenuation  mechanisms  must  also  be  considered.  These  losses 
are  caused  by  viscosity,  thermal  conductivity,  mass  diffusion,  and  vibra¬ 
tional  and  rotational  relaxation.  In  addition,  there  are  losses  associated 
with  sound  propagation  inside  a  tube.  Attenuation  resulting  from  viscosity 
and  thermal  conductivity  is  termed  classical  attenuation.  The  classical 
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attenuation  coefficient  is  given  by 

aCL  »  21r2/pa2(|n+  ^^)f2  (56) 

where  ac^  is  the  exponential  decay  constant  for  sound  pressure  amplitude, 
i.e.  Aexp(-ax)  where  x  is  the  propagation  distance.  Beginning  with  the 
conservation  equations  for  mass,  energy,  and  momentum,  Kirchhoff21  obtained 
an  algebraic  equation  for  the  propagation  constant  for  radially  summetric 
waves.  Kirchhoff  assumed  that  particle  velocity  and  sound  temperature  were 
zero  at  the  tube  wall.  Henry22  pointed  out  that  these  conditions  might  be 
inadequate,  owing  to  the  temperature  jump  and  slip  velocity  at  the  tube 
wall.  Shields23  has  modified  Kirchhoff 's  equations  to  include  temperature 
jump  and  slip  velocity.  The  corrected  equations  were  solved  numerically 
using  a  computer  program  developed  by  Shields  which  is  described  in  Refer¬ 
ence  24.  The  physical  constants  for  our  waveguide  were  obtained  from 
Shields25.  Reflection  coefficients  at  the  ends  of  the  tube  were  needed  to 
correct  for  losses  upon  reflection.  The  coefficients  were  measured  at  room 
temperatures.  The  temperature  dependence  of  reflection  coefficients  is 
weak21**26  and  therefore  was  not  considered;  endplate  temperature  varied 
only  slightly  during  our  experiments. 

When  a  density  gradient  exists  in  a  mixture  of  gases,  there  will  be 
a  greater  flow  of  light  molecules,  relative  to  their  concentration,  than  of 
heavy  molecules.  This  results  because,  for  a  given  temperature,  lighter 
molecules  are  traveling  faster  than  heavy  molecules.  The  results  is  a 
periodic  unmixing  of  the  gas,  offset  by  irreversible  diffusion,  which 
results  in  an  energy  loss  of  the  sound  wave.  This  effect  is  additive  to 


vV 
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that  of  classical  absorption  and  has  an  attenuation  coefficient  given  by 

s 

adiff  *  (TT2yf2/a3p2)  £  xjSit(mj-mk)x(mjDij-micDik)  (57) 

i.j.k 

where  D^j  is  the  multicomponent  diffusion  coefficient  for  particle  "i"  and 
"j",  m±  the  mass  of  the  i^1  particle.  For  the  gas  mixtures  used  in  this 
study,  attenuation  resulting  from  diffusion  was  much  smaller  than  those 

caused  by  tube  losses  or  vibrational  relaxation.  The  computer  programs 
which  calculated  attenuation  did  not  consider  attenuation  resulting  from 
diffusion. 

When  sound  propagates  through  polyatomic  gases  the  compressions  and 
rarefactions  involve  a  redistribution  of  energy  between  the  various  modes 
of  excitation;  this  includes  vibrational  and  rotational  modes.  At  suffi¬ 
ciently  high  frequencies,  the  period  of  fluctuation  becomes  shorter  than 
the  time  required  for  the  redistribution  of  energy  into  the  various  modes- 
When  this  occurs  the  compressions  and  rarefactions  are  no  longer  reversible 
and  attenuation  results.  Since  rotational  energy  levels  are  closely  spaced 
near  room  temperature  for  CI2 ,  H2 ,  and  SF6 ,  a  single  relaxation  time  can  be 
used27  to  describe  rotational  relaxation.  The  resulting  absorption  coeffi¬ 
cient,  for  the  frequencies  we  used,  is27 
s 

arot  ■  {[Nx{irf2(di/2)R(Y-l)]/(fr,rotc?a)i  (58) 

i 

where  d^  is  the  number  of  rotational  degrees  of  freedom,  i.e.  2  for  H2 , 

CI2,  3  for  SF6;  xi  the  fractional  number  of  molecules  of  species  "i";  N 
is  the  total  number  of  molecules;  and  ff,rot  the  frequency  associated  with 
relaxation  of  rotational  modes.  Since  arot  and  have  the  same  frequency 
dependence,  rotational  relaxation  can  be  included  by  modifying  the  classi¬ 
cal  attenuation  coefficient26. 
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Two  vibrational  attenuation  mechanisms  were  considered  in  our  calcu¬ 
lations:  collisions  which  relax/excite  the  first  vibrational  level  in  Cl2 ; 
and  collisions  which  relax/excite  vibrational  levels  of  SF6.  H2  and  upper 
Cl2  vibrational  states  are  not  considered  because  at  lower  temperatures 
they  are  rarely  excited.  The  first  mechanism  is  especially  simple  because 
it  involves  only  one  vibration-translation  (v-t)  process  and  can  be  consid¬ 
ered  a  pseudo-first-order  chemical  reaction5 

Cl2 (v*0)  t  Cl2 (v*l)  (59) 

In  general,  collisions  which  exchange  multiple  quanta  of  energy,  in  v-t  and 
v-v  processes,  must  be  considered  to  predict  vibrational  attenuation.  For 
example,  if  A  and  B  collide,  A  may  loose  two  quanta  of  vibrational  energy 
and  B  receive  one  quantum  of  vibrational  energy  and  some  translational 
energy.  This  coupling  between  v-t  and  v-v  processes  complicates  the  prob¬ 
lem28*29  greatly;  however,  SF6  has  properties  which  simplify  the  problem. 
SF6  has  many  vibrational  levels  which  are  excited  even  at  low  temperatures 
(~300K)  but,  the  lowest  level  deexcites  much  slower  than  the  others.  This 
lowest  level  gating  makes  it  possible  to  characterize  SF6  vibrational 
relaxation  with  one  relaxation  time29,  Tvt  •  For  frequencies  such  that 
UTvt  K<  1  »  which  includes  the  range  we  used , 

<60> 

where  CQ  is  the  constant  volume  specific  heat  for  low  frequencies,  and 
Cv  is  the  contribution  to  Cc  resulting  from  vibrational  relaxation. 

Eq.  (60)  is  multiplied  and  divided  by  pressure,  p,  because  (pxvt)  is 


•<? 


pressure  independent  for  binary-collision  induced  relaxation.  The  vibra¬ 
tional  relaxation  time,  tvt,  for  a  mixture  of  CI2-H2-SF6  is  given  by 

1  _  x-H2  ,  x'Clj  (  x^SFs 

Tvt  Tvt,H2  Tvt,Cl2  Tvt,SF6 

where  Tvt,SF6is  the  vibrational  relaxation  time  for  pure  SF6 ;  Tvt,j  is  the 
vibrational  relaxation  time  related  to  SF6  in  mixtures  of  "j".  Equation 
(61)  can  be  used  to  determine  the  constants  Tvt,SF6*  Tvt  Cl2’  an<*  Tvt  H2 
determining  Tvt  for  a  number  of  known  concentrations,  x'g^,  x^ci2»  ant* 
x"sf6.  TheTvt15  are  found  by  measuring  atotal  for  various  pressures, 
frequencies,  and  concentrations  of  H2  and  Cl2  in  SF6;  substracting  other 
attenuation  terms  to  get  ei^  for  each  set  of  parameters;  and  solving  eq. 
(60)  for  xvt. 

Reaction  Details 

Two  SACERs  were  constructed  for  this  investigation;  one  operated  in  a 
pulsed  mode,  the  other  in  a  cw  mode.  Both  systems  are  characterized  chem¬ 


ically  by  the  following  H2-C12  reactions3’11* 

Cl2  +  hv  ■+■  2C1  (photo  initiated  dissociation)  (62) 
Cl  +  H2  t  HC1  +  H  (63) 
H  +  CI2  *  HC1  +  Cl  (64) 
M  +  2CL  *  M  +  Cl2  (65) 
Cl2(v-0)  t  Cl2(v-1)  (66) 


Where  M  represents  SF6  for  the  cw  SACER,  Ar  for  the  pulsed  SACER.  The  last 
"reaction"  is  a  pseudo-first-order  reaction  which  allows  for  CI2  vibrational 
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CHAPTER  III 

APPARATUS 

Two  systems  were  constructed  for  this  investigation.  The  first  was 
used  for  pulsed  SACER  experiments,  the  second  for  cw-SACER  and  attenua¬ 
tion  experiments.  A  brief  description  of  the  systems  will  be  given, 
followed  by  a  more  complete  discussion  of  individual  components. 

Fig.  (1)  illustrates  the  acoustic  and  lighting  equipment  used  for 
the  cw  experiments.  The  heart  of  the  system  was  a  cylindrical  waveguide 
(lm  x  3cm)  with  transducers  at  each  end.  The  tube  was  made  of  Vycor 
which  is  UV-transparent .  Three  UV-fluorescent  lamps  surrounded  the 
tube  and  provided  UV-radiation  needed  to  dissociate  CI2 .  Aluminum  slats 
supported  the  lamps  and  helped  insure  homogenous  illumination  of  the 
Vycor  tube.  The  acoustic  waveguide  and  lamps  were  enclosed  in  an  alumi¬ 
num  tube  to  Increase  UV-intensity  inside  (aluminum  is  an  excellent  UV- 
ref lector) .  The  entire  system  was  covered  with  a  light-tight  cloth  to 
protect  the  lab  from  UV-radiation  and  isolate  the  system  from  UV-radia¬ 
tion  emitted  by  fluorescent  lamps  in  the  lab.  A  similar  system  was  used 
for  pulsed  experiments  (see  fig.  (2)).  Intense  flash  lamps  were  used  to 
dissociate  CI2  at  the  start  of  each  experiment.  Unilluminated  regions 
at  each  end  of  the  tube  generated  acoustic  pulses  which  propagated  down 
the  tube.  Two  microphones  were  used  to  monitor  the  amplitude  of  these 
pulses.  The  first  microphone  was  mounted  flush  in  one  end  of  the  tube, 
the  second  microphone  was  mounted  near  the  other  end  of  the  tube.  A 


cw  System 


Fig.  2 

Pulsed  SACER 
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movable  sending  transducer  was  mounted  opposite  the  first  microphone, 
but  was  used  only  as  a  reflector  because  the  intense  acoustic  pulses, 
generated  by  the  unilluminated  regions,  masked  tone  bursts  generated  by 
the  transducer.  A  LSI-11  based  minicomputer  (MINC  system  11/23)  was 
used  for  cw-SACER  experiments  to  receive  data,  control  experiment,  and 
graph  and  analyze  data.  The  computer  was  equipped  with  an  A/D  converter 
and  digital  out  (DO)  module.  The  A/D  converter  was  used  to  digitize  the 
signal  from  the  receiver.  The  DO  module  generated  digital  pulses  which 
controlled  the  sending  transducer  and  equipment  associated  with  fluores¬ 
cent  lamps.  Attenuation  measurements  were  made  with  the  system  used  for 
cw  experiments.  The  lamps  were  not  needed  for  attenuation  measurements. 
All  experiments  used  research  grade  gases  (Matheson  Gas,£  99.96 %  pure) 
with  the  exception  of  H2  used  during  attenuation  measurements.  Low 
impurity  levels  could  be  tolerated  in  H2-SF6  attenuation  measurements 
because  H2  is  generally  more  effective  at  relaxing  vibrational  modes  in 
SF6  than  impurities. 

Two  solid  dielectric  capacitance  transducers  were  used  in  the  cw 
system  (see  fig.  (3)).  Both  were  constructed  in  our  laboratory  for  this 
experiment.  Each  had  a  high  voltage  plate  covered  with  0.0127  mm  Kapton 
which  had  an  aluminum  coating  on  the  side  away  from  the  high  voltage  back 
plate.  The  aluminum  coating  and  back  plate  formed  a  capacitor  whose 
separation  distance  could  be  varied  with  acoustical  pressure  or  applied 
voltage.  One  transducer  was  operated  in  transmitting  mode,  the  other  in 


a  receiving  mode. 


TABLE  (1) 


Equipment  List 

1.  General  Radio  Company,  Preamplifier  Type  1560-p42 

2.  General  Electric,  Black  Light  no.  F40BL 

3.  Hewlett-Packard,  3312A  Function  Generator  for  cw-SACER  experiments; 
Hewlett-Packard,  200  CDR  Wide  Range  Oscillator  for  attenuation 
experiments 

4.  Dynaco,  Mark  VI  Amplifier 

5.  Electro  Instruments,  Inc.,  Precision  Power  Source  no.  6303PS 

6.  General  Radio  Company,  1560-P62  Power  Supply 

7.  Rockland,  Dual  Hi/Lo  Filter  no.  452 

8.  Eg&G  Parc,  113  Pre-Amp 

9.  Each  lamp  had  its  own  Regulated  Power  Supply: 

Philbrick  Researchers,  R-600 
Philbrick  Researchers,  R-300 
Lambda,  C-881M 

10.  Sheldahl,  0.0127  mm  Kapton  with  aluminizing 

11.  Fluke,  1911-A  Multi-Counter 

12.  Eico,  1064  Power  Supply 

13.  Digital,  minicomputer  no.  MNC 11-FA,  version  2.0,  system  11/23, 

LSI- 11  based 

14.  Albia  Electronics,  DM-6  Variable  Power  Supply 

15.  Xeon  Corporation,  Model  A  Power  Supply 

16.  Xeon  Corporation,  Model  C  rf-Trigger 

17.  Xeon  Corporation,  1  meter  Micropulse  Flash  Lamp 

18.  Bruel  &  Kjaer,  Condenser  Microphone  Cartridge  Type  4149 
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For  the  transmitter,  a  force  was  applied  to  the  aluminum  coating  by 
varying  the  potential  between  it  and  the  high  voltage  plate.  This  force 
caused  the  Kapton  to  oscillate  sending  an  acoustic  signal  down  the  tube. 
Two  voltages  were  applied  to  the  transmitter,  a  200  volt  bias  which 
insured  that  the  foil  and  back  plate  did  not  separate,  and  a  ±  20  volt 
tone  burst  which  caused  the  oscillations.  The  two  voltages  were  combined 
with  a  voltage  adder  built  in  this  lab  (see  appendix  B) .  The  sender  was 
movable  so  the  distance  between  transmitter  and  receiver  could  be  changed. 

Acoustical  pressure  variations  on  the  receiving  transducer  caused 
the  potential  between  aluminum  coating  and  backplate  to  vary.  A  GenRad 
preamp/power  supply  combination  was  used  with  the  receiver  to  apply  a 
200  volt  bias  and  detect  the  ac  voltage  caused  by  acoustic  signals.  The 
receiver  sensitivity  was  controlled  by  decreasing  the  active  area  of  the 
high  voltage  plate  (see  fig.  (4)).  Fig.  (5)  illustrates  additional 
signal  conditioning  electronics  associated  with  both  transducers. 

The  pulsed  system  was  equipped  with  two  1.27  cm,  B&K  condenser 
microphones  (see  fig.  (2)).  One  microphone  was  mounted  flush  in  an 
endplate,  the  other  was  mounted  near  the  other  endplate  with  its  active 
surface  parallel  to  the  tube's  axis.  Both  were  powered  by  GenRad  preamp/ 
power  supply  combination  like  the  one  described  for  the  cw  system.  A 
sending  transducer  was  mounted  in  the  other  endplate;  the  sound  it  gener¬ 
ated  was  not  intense  enough  to  be  used  for  the  pulsed  experiments. 

Several  constraints  were  placed  on  the  UV-radiation  sources  used  for 
cw  operation.  They  had  to  start  quickly,  relative  to  the  relaxation 
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times  of  the  chemical  reactions,  and  not  cause  large  electrical  disturb¬ 
ances  which  would  be  picked  up  by  the  receiver  preamp.  To  this  end,  a 
complicated  procedure  was  developed  to  start  the  UV-f lourescent  lamps 
(see  fig.  (6)): 

1.  300  volts  dc  was  applied  across  the  lamps 

2.  12  volts  dc  was  supplied  to  glow  one  filament  for 
1/2  sec. 

3.  the  lamps  started  when  the  wire,  wrapped  around  the 
lamps,  and  the  slats  were  shorted  to  the  300  volt 
terminal;  the  wire  and  slats  were  "floating"  before 
the  shorting  process. 

After  starting,  the  lamps  drew  one  amp  each.  Due  to  the  large  power 
requirement,  each  lamp  was  powered  separately.  Timing  for  lamp  starting 
and  transducer  activation  was  provided  by  the  minicomputer.  The  UV-flash 
lamps  used  for  the  pulsed  experiments  had  to  dissociate  the  Cl2  quickly, 
compared  to  the  time  scale  of  the  chemical  reactions,  at  the  beginning 
of  each  experiment.  The  four  cylindrical  flash  lamps  (lm  x  1cm)  were 
powered  by  a  bank  of  capacitors.  They  were  triggered  with  a  rf-oscil- 
lator  (see  fig.  (7)).  The  lamps  were  flashed  at  10  kvolts  with  a  30 
kvolt  trigger.  Pulse  duration  was  estimated  by  the  supplier  (Xenon,  Inc.) 
to  be  10  psec.  The  intensity  of  the  UV-flash  could  be  changed  by  varying 
the  number  of  capacitors  in  the  power  supply. 

As  mentioned  earlier,  the  cw-SACER  experiments  were  controlled  by  a 
MINC-11/23  computer  system  manufactured  by  Digital  Equipment  Corporation. 
The  system  included  lab  modules  which  digitized  data  and  provided  digital 
signals  to  activate  lab  equipment.  Fig.  (8)  illustrates  the  computer  and 
related  equipment.  The  digital  output  (DO)  module  sent  digital  pulses, 
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Electronics  for  UV-Lamps  with  cw-SACER 


Electronics  for  UV-Flash  Lamps  Used  in  Pulsed  System 


Minicomputer  System  Used  to  Control  cw-SACER 
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of  variable  length  and  time,  to  the  electronic  buffer  equipment.  The 
buffer  equipment  insured  that  DO  signals  were  impedance-matched  with 
later  equipment  and  that  electronic  noise,  caused  by  the  SACER,  did  not 
get  sent  back  to  the  computer.  Three  DO  channels  were  used  to  control 
the  function  generator,  current  to  lamp  filaments,  and  voltage  applied 
to  the  wire  wrapped  around  the  lamps.  The  A/D  converter  was  used  to 
digitize  the  amplified  signal  from  the  receiving  transducer. 

During  cw-SACER  experiments  the  amplitude  of  a  tone  burst  was 
monitored  as  it  propagated  through  a  reacting  gas  mixture.  UV-lamps 
continuously  dissociated  Cl 2  allowing  the  reaction  to  proceed  slowly 
until  all  the  H2  was  used  (approximately  100  msec) .  The  experiments 
were  conducted  as  follows: 


1.  Using  the  law  of  partial  pressures,  a  stainless 
steel  storage  tank  was  filled  with  enough  Cl  2, 

H2,  and  SF6,  in  the  desired  concentrations,  to 
conduct  several  experiments.  The  gases  were  then 
allowed  to  mix  one  day  for  each  atmosphere  of  pres¬ 
sure  in  the  tank. 

2.  After  mixing,  but  before  the  SACER  tube  was  filled 
with  gas,  the  lamps  were  turned  on  and  off  once  to 
insure  proper  starting  during  later  operation. 

3.  300  volts  was  applied  across  the  lamps  without 
starting  them. 

4.  The  SACER  tube  was  filled  with  mixed  gas. 

5.  The  computer  then  controlled  the  SACER  in  this  order: 
-glow  lamp  filaments  1/2  sec  to  insure  proper  lamp 

starting 

-then  start  digitizing  data  from  receiving 
transducer  and  start  lamps  by  changing  the  poten- 
of  the  wire  coiled  around  the  lamps. 

-activate  sending  transducer  to  send  4  msec  tone 
burst  at  predetermined  time  after  lamp  start 
-stop  digitizing  130  msec  after  lamp  start. 

6.  After  each  run  the  lamps  were  turned  off  and  the 
gases  pumped  out  of  the  SACER  tube. 
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7.  Seeps  2-6  were  repeated  at  approximately  half 
hour  intervals  until  the  storage  tank  was  de- 
pleated. 

During  each  pulsed  experiment,  UV-flash  lamps  dissociated  Cl2  in  the 
SACER  tube.  This  began  the  reaction  and  caused  the  pressure  to  rise 
rapidly  everywhere  inside  the  tube  except  for  small  unilluminated  re- 
tions  at  each  end  of  the  tube.  The  increasing  pressure  caused  gas  to 
expand  into  the  unilluminated  regions.  This  expansion  led  to  pulses 
which  propagated  down  the  tube.  Microphones  monitored  the  amplitude 
of  the  pulses  as  they  propagated  inside  the  tube.  Because  all  the  Cl2 
was  immediately  dissociated,  the  reactions  only  lasted  approximately 
10  msec.  Pulsed  SACER  experiments  were  conducted  as  follows: 

1.  The  SACER  tube  was  filled  with  the  proper  mixture 
of  Cl2 ,  H2 ,  and  Ar.  The  gases  were  then  allowed 
to  mix  thoroughly  (at  least  four  hours) . 

2.  The  capacitors  were  charged  and  then  used  to  flash 
the  lamps.  The  microphone  outputs  were  stored  on 
a  fm- recorder  and  later  analyzed. 

In  determining  the  vibrational  relaxation  time  for  gas  mixtures 
used  during  cw-SACER  experiments,  a  set  of  attenuation  experiments  were 
conducted.  During  attenuation  experiments,  the  amplitude  of  the  first 
reflection  was  measured  for  different  transducer  separations.  Experi¬ 
ments  were  conducted  at  various  pressures,  frequencies,  and  Cl2  and  H2 
concentrations  in  SF6 .  The  concentrations  were  10,  20,  and  50  percent 
H2 ,  10  and  20  percent  Cl2 ,  and  pure  SF6 ;  the  frequencies  ranged  from 
1.5-50kHz;  and  the  pressures  ranged  from  25-400  torr.  These  steps  were 
followed  for  each  attenuation  experiment: 


A  stainless  steel  storage  tank  was  filled  with 
enough  SF6  and  CI2  or  H2 ,  in  the  desired  partial 
pressures,  to  fill  the  tube  several  times. 

After  mixing  the  gases  for  24  hours,  the  SACER 
tube  was  filled  to  the  highest  pressure  used  in 
attenuation  experiments. 

At  the  largest  transducer  separation  distance, 
the  pressure  amplitude  of  the  first  received  tone 
burst  was  measured  for  each  frequency. 

The  distance  between  transducers  was  reduced  and 
step  3  repeated. 

Step  4  was  repeated  approximately  twelve  times. 
The  pressure  was  lowered  and  steps  3-5  repeated. 
Because  a  small  amount  of  leakage  occurred  around 
the  movable  transducer's  mount,  the  system  was 
periodically  evacuated  and  refilled  to  ensure  the 
gas  mixture  was  not  contaminated. 


CHAPTER  IV 


EXPERIMENTAL  RESULTS 

Pulsed  SACER 

Pulsed  SACER  experiments  were  begun  by  generating  Cl  atoms  with 
a  very  short  UV  pulse,  then  allowing  reactions  (63-66)  to  go  to  comple¬ 
tion  modulated  by  an  acoustic  signal.  The  acoustic  signal  resulted 
from  unilluminated  regions  in  the  tube  and  was  monitored  with  two 
microphones.  Electronic  filters  were  used  to  examine  individual  fre¬ 
quency  bands  found  in  the  acoustic  signal.  A  frequency  band  centered 
about  2kHz  was  examined  for  amplification. 

Experimentally,  the  amounts  of  Cl2 ,  H2 ,  and  Ar  initially  in  the 
system  could  be  controlled.  Also,  varying  the  flash  intensity  gave  some 
control  over  the  initial  concentration  of  atomic  chlorine.  One  of  the 
first  tasks  was  to  determine  the  optimum  conditions  for  observing 
sound  amplification.  It  was  necessary  to  consider  the  chemical  reactions 
(63-65)  and  the  attenuation  mechanisms  present.  A  computer  program  writ¬ 
ten  by  Gilbert  et  al.5  was  used  to  include  the  chemical  reactions;  a 
program  developed  by  Shields23  was  modified  to  compute  tube  and  classical 
attenuation  as  a  function  of  frequency,  gas  composition,  temperature,  and 
pressure.  These  two  programs  were  used  to  solve  for  signal  strength  as 
a  function  of  time  for  various  reacting  gas  mixtures. 

Ideally,  we  would  like  to  have  computed  amplication  for  many  sets 
of  initial  conditions;  however,  since  each  set  required  about  30  CPU 
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minutes  (Dec  system  1077)  only  60  sets  were  considered.  From  the  compu¬ 
ter  outputs  we  were  able  to  develop  a  physical  feel  for  the  competing 
mechanisms.  For  example,  Table  (2)  illustrates  the  effect  of  varying 
H2  concentration.  For  these  calculations,  it  was  assumed  that  0.0005 
atm  of  CI2  were  dissociated  by  the  UV  flash.  At  low  H2  concentrations, 
the  primary  mechanism  for  amplif iciation  is  recombination  of  CI2.  In 
the  absence  of  tube  and  classical  losses,  after  3msec,  this  mechanism 
would  result  in  a  gain  of  1.7.  Attenuation  decreases  this  to  1.5  and 
0.75  at  lm  and  2m  respectively.  As  the  H2  concentration  is  increased, 
reaction  (64)  becomes  more  important  and  the  amplification  increases  up 
to  about  8  after  lm  because  reaction  (64)  is  highly  exothermic.  Greater 
H2  concentration  does  little  to  increase  amplification. 

Table  (2)  also  illustrates  that,  for  the  case  of  varying  H2  concen¬ 
tration,  large  amplification  is  accompanied  by  large  overpressure.  This 
trend  occurred  repeatedly  in  our  computer  simulations.  The  microphones 
used  in  the  pulsed  system  could  only  withstand  an  overpressure  of  0.03 
atm  which  limited  allowable  initial  conditions  and  expected  amplification. 

In  pulsed  SACGR  experiments,  after  the  lamps  dissociate  Cl2,  the 
pressure  in  the  illuminated  region  increases  due  to  the  dissociation 
and  the  subsequent  heating  from  exothermic  reactions.  The  increased 
pressure  in  the  illuminated  region  expands  into  the  unilluminated  regions 
at  each  end  of  the  tube.  When  the  expansion  reaches  the  microphones,  the 
pressure  rapidly  increases.  The  positive  pressure  pulse  is  followed  by  a 
decrease  in  pressure.  The  lower  pressure  duration  is  about  Xn/c  where 
c  is  the  speed  of  sound  and  Xn  is  the  length  of  the  unilluminated  section. 


TABLE  (2) 


Pulsed  SACER  Results 


Total  Pressure  *  0.025  atm 


Concentrations  In  atm  Overpressure  (atm) 


Amplitude 
(initially  1) 
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This  pattern  is  repeated  as  the  pressure  pulses  bounce  back  and  forth 
in  the  tube.  An  example  of  this  type  of  signal  is  illustrated  in  figure 
(9).  By  filtering  the  signal,  the  amplitude  of  one  band  of  frequencies 
can  be  observed.  When  this  is  done,  a  series  of  decaying  pulses  is 
observed;  each  pulse  being  similar  to  a  sine  wave  with  a  gaussian 
envelope.  It  is  the  change  in  amplitude  of  these  pulses  that  is  examined 
for  amplification. 

The  system  had  two  microphones  which  detected  the  pressure  changes. 
One  microphone  was  mounted  flush  in  the  end  opposite  the  sending  trans¬ 
ducer,  the  other  was  mounted,  with  its  face  parallel  to  the  axis  of  the 
sound  tube,  at  the  same  end  as  the  sending  transducer.  These  microphones 
will  be  referred  to  as  mike  1  and  mike  2  respectively. 

For  the  optimum  mixture  listed  above,  the  first  burst  arrives  at 
mike  1  about  2msec  after  the  flash.  The  second  burst  has  reflected  off 
the  opposite  end,  traveled  back  down  the  tube  and  arrives  about  7msec 
after  the  flash.  The  third  burst  is  the  same  as  the  first  after  traveling 
to  mike  2,  reflecting  off  the  sending  transducer  and  returning.  Subse¬ 
quent  pulses  correspond  to  reflections  back  and  forth  in  the  tube. 

For  the  purpose  of  the  following  arguments,  a  numbering  system  has 
been  constructed  to  keep  track  of  these  bursts.  The  first  burst  received 
at  mike  1  is  denoted  Bj,  the  second  B^,  the  third  B*,  etc.  At  mike  2,  the 
same  bursts  are  referred  to  as  B§,  B*,  B|,  etc.  In  this  way,  an  odd  sub¬ 
script  refers  to  a  burst  which  arrived  first  at  mike  1.  This  means  that 
Bj[  is  the  same  burst  as  B^  but  B^  has  traversed  an  additional  lm  down  the 
tube  to  mike  2.  The  amplitude  of  B^  will  be  denoted  by  A*,  etc. 


*  *  **  »  "  «  ‘  m  *  »  *  .  * 


Before  discussing  specific  results,  it  should  be  noted  that  the 

ratio,  Aj/Aj  has  no  particular  significance  by  itself  since  this  ratio  is 

a  function  of  the  unilluminated  volumes,  on  each  end,  which  (for  the  mike 

2  end)  is  variable.  Further,  since  mikes  1  and  2  are  mounted  differently, 

the  ratio  Aq/aJ  has  no  significance.  However,  by  comparing  several  experi¬ 
ments  these  ratios  can  provide  information. 

Figure  (10)  shows  the  signal  strength  at  mike  1,  on  a  log  scale,  ver¬ 
sus  propagation  distance  for  Cl2-Ar  mixture.  For  a  gas  with  no  amplifica¬ 
tion,  the  points  should  lie  along  a  straight  line  since  the  signal  strength 
A  is  equal  to  Ao  exp(ax-ng)  where  a  is  the  absorption  coefficient  for  the 
gas,  B  is  the  reflection  coefficient,  and  n  the  number  of  reflections.  As 
shown  in  figure  (10),  the  points  lie  along  a  straight  line  indicating  there 
is  little  if  any  amplification  from  the  Cl2-Ar  mixture,  as  expected.  Also 
note  that  although  the  points  for  the  signal  originating  near  mike  2  are 
offset  in  absolute  magnitude,  the  decay  rate  is  the  same  (to  within  the 
measurement  accuracy).  Figure  (10)  indicates  that  the  2kHz  component  of 
the  pressure  pulse  originating  near  mike  2  begins  at  a  level  1.15  greater 
than  the  pulse  originating  near  mike  1. 

Figure  (11)  illustrates  the  results  of  adding  H2  to  the  Cl2-Ar 
mixture  in  the  ratio,  0.077  atm  Ar:  0.019  atm  CI2:  0.003  atm  H2.  The 
signal  strengths  for  both  microphones  (obtained  from  oscilloscope  traces 
of  the  filtered  signals)  are  plotted  on  a  log  scale.  At  both  microphone 
positions,  the  signal  decreases  logrithmically  with  distances,  the 
recorded  amplitude  is  well  below  the  straight  line  extrapolation.  This 
suggests  that  there  is  a  time  over  which  the  signal  is  amplified.  By 


Figure  10 

Signal  Strength  Versus  Propagation  Distance  for  2kHz  Signal 
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extrapolating  the  line  back  to  1/2  m,  we  determined  the  effective  gain 
which  took  place  in  the  time  necessary  for  the  signal  to  travel  this 
3/4  m;  the  gain  was  1.8  (two  other  sets  of  data  give  gains  of  1.8  and 
1.5).  This  means  that  for  a  propagation  distance  greater  than  3/4  m, 
the  signal  amplitude  is  1.8  times  what  it  would  have  been  without  ampli¬ 
fication.  A  gain  of  4  was  predicted  for  the  mixtures  used. 

There  are  additional  features  of  this  data  which  serve  to  convince 
us  that  amplification  has  been  observed  in  the  pulsed  system.  Comparing 
figures  (10)  and  (11) ,  note  that  for  mike  1  with  H2 ,  B*  is  larger  in 
magnitude  than  Bi;  in  the  absence  of  H2 ,  b£  is  larger.  The  geometry  was 
the  same  so  the  ratio  of  A*  to  A2  should  be  the  same  unless  B*  has 
experienced  more  amplification  in  the  1/4  m  greater  distance  it  traveled 
compared  to  Bj. 

Both  these  observations  seem  convincing,  but,  unfortunately  the 
results  are  not  conclusive.  The  large  unfiltered  pulse  amplitude  raises 
the  possibility  of  non-linear  effects  in  the  gas.  One  would  generally 
expect  a  non-linear  wave  to  decay  more  rapidly  leading  to  a  curve  with 
decreasing  steepness,  opposite  to  that  seen  in  figure  (11).  But  the 
fundamental  of  the  pulse  is  at  a  much  lower  frequency  so  there  could  be  a 
build  up  at  higher  frequencies  due  to  steepening  of  pulse  wave  forms. 

The  appearance  of  steepening  is  not  obvious  from  the  data,  yet  we  are  not 
prepared  to  rule  out  this  possibility  entirely. 

The  magnitude  of  the  observed  gains  are  about  1/2  that  predicted  by 
Gilbert’s  computer  program.  There  are  several  possible  explanations  for 
this  discrepancy.  The  computer  program  did  not  include  wall  collisions 
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or  expansion  into  the  unilluminated  regions  (which  cools  the  gas). 

Wall  collisions  or  expansion  could  slow  the  reaction  scheme.  The  effect 
of  wall  collisions  should  be  small  and  could  be  explored  by  reproducing 
the  experiment  at  lower  pressures  where  wall  collisions  are  more  impor¬ 
tant.  Cooling  effects  caused  by  the  unilluminated  regions  were  overcome 
in  the  cw-SACER  experiments. 

Another  possible  cause  for  the  difference  between  observed  and 
predicted  amplification  rates  is  the  uncertainty  in  the  UV  energy  sup¬ 
plied  by  the  flash  lamps.  The  energy  from  the  flash  lamps  was  estimated 
by  flashing  (with  a  given  flash  intensity)  with  different  Cl2  pressures 
in  the  tube,  observing  the  total  overpressure,  then  comparing  this  to 
the  theory  with  the  original  Cl  produced  used  as  an  adjustable  variable 
to  obtain  agreement  between  theory  and  prediction.  An  example  is  given 
in  figure  (12).  This  technique  neglects  the  unilluminated  regions  and, 
hence,  underestimates  the  Cl  concentration  in  the  illuminated  region.  But 
that  underestimate  should  not  be  greater  than  10  percent  (the  volume 
ratio)  which  would  not  entirely  explain  the  difference  between  observed 
and  predicted  gain. 

cw-SACER 

During  cw-SACER  experiments,  the  UV-flourescent  lamps  continuously 
dissociated  Cl2,  providing  atomic  Cl  for  reactions  62  through  66.  For  all 
but  one  experiment,  a  tone  burst  was  sent  down  the  tube  108msec  after  the 
lamps  were  started;  the  receiving  transducer  monitored  the  tone  burst  as 


it  propagated  back  and  forth  inside  the  tube.  Tone  bursts  were  also  sent 
before  and  after  each  experiment  as  a  control. 

Experimentally,  the  amounts  of  Cl2,  H2,  and  SF6,  initially  in  the 
system,  could  be  controlled.  Unlike  the  pulsed  system,  lamp  intensity 
was  not  variable.  Initial  gas  concentrations  similar  to  those  used  by 
Toong7  were  selected  for  the  cw-SACER  experiments.  Three  C12:H2:SF6 
ratios,  at  a  pressure  of  1/3  atm,  were  selected,  0.5:2:17.5,  1:2:17, 
and  2:2:16.  For  each  ratio,  experiments  were  conducted  for  four  acoustic 
frequencies,  1.0,  2.5,  4.0,  and  6.5kHz. 

For  the  Initial  conditions  selected,  the  reactions  proceeded  slowly 
and  gave  predicted  amplification  which  varied  only  a  small  amount  in  the 
time  it  took  for  a  tone  burst  to  decay  to  below  the  noise  level.  There¬ 
fore  the  tone  bursts,  sent  while  the  reactions  were  proceeding,  decayed 
exponentially  with  distance,  but  with  a  smaller  attenuation  coefficient 
than  tone  busts  sent  when  the  reactions  were  not  progressing.  The 
filtered  signal  from  the  receiving  transducer  was  digitized  and  later 
used  to  obtain  plots  of  tone  burst  amplitude  versus  distance  and  time 
(see  fig.  (13)). 

Table  (3)  and  fig.  (14)  show  the  data  from  the  cw-SACER  experiments 
just  described.  The  plots  of  tone  burst  amplitude  versus  distance  and 
time  were  used  to  obtain  attenuation  coefficients  and  velocities.  In 
table  (3),  the  quantities  listed  in  the  "before"  columns  are  averages  of 
measurements  taken  before  and  after  the  reactions.  After  the  lamps  were 
turned  off  and  the  mixture  was  allowed  to  cool,  the  attenuation  coefficient 
and  sound  velocity  were  the  same  as  before  the  reactions,  within 


TABLE  (3) 


cw  SACER  Data 


Frequency 

and  Attenuation 

Mixture  (dB/cm)  Aa 


Velocity 

(m/sec) 


H2-CI2-SF6  Before  During  (dB/cm)  Before  Durin; 


2-2-16 

1.0  kHz 

.0106 

1 

.0019 

■ 

1 

2.5 

.0170 

.0014 

■PPJ 

4.0 

.0240 

.0024 

156 

■H 

6.5 

.0371 

.0032 

155 

165 

1-2-17 

1.0  kHz 

.0108 

.0098 

.0010 

152 

155 

2.5 

.0177 

.0162 

.0015 

154 

156 

4.0 

.0237 

.0212 

.0025 

151 

156 

6.5 

.0357 

.0329 

.0028 

152 

156 

0.5-2-17.5 
1.0  kHz 

.0107 

.0102 

.0005 

150 

151 

2.5 

.0174 

.0168 

.0006 

149 

152 

4.0 

.0238 

.0223 

.0006 

150 

152 

6.5 

.0356 

.0346 

.0010 

149 

153 

experimental  error.  Aa  is  the  change  in  attenuation  while  the  reactions 
were  taking  place  and  is  obtained  by  subtracting  the  "during"  attenuation 
coefficient  from  the  "before"  attenuation  coefficient.  AT  is  the  dif¬ 
ference  between  the  temperature  when  the  reactions  are  progressing 
(determined  from  sound  speed  when  the  tone  burst  is  decaying)  and  room 
temperature.  AT  is  computed  from  the  velocity  measurements  by  using  the 
following  eq.  derived  from  the  equation  for  adiabatic  sound  velocity. 

Wing  -  W.  uSfflEiSS)  (67) 

vBefore 

The  data  from  the  2:2:16  mixture  was  compared  with  Gilbert's  computer 
program6;  tube  losses  and  classical  attenuation  were  included  in  this 
calculation  by  using  a  slightly  modified  computer  program  originally 
developed  by  Shields23.  The  frequency  factor,  for  reaction  (62),  used 
in  the  computer  predictions  was  that  value  which  correctly  predicted  the 
temperature  change  during  the  reactions.  This  frequency  factor  (see  first 
assumption  under  eq.  (49))  is  related  to  the  intensity  of  the  fluorescent 
lamps.  Table  (4)  and  fig.  (14)  present  the  comparison  between  the  pre¬ 
dicted  and  observed  attenuation  coefficients.  The  observed  attenuation 
coefficients  are  in  fair  agreement  with  the  computer  predictions.  Also, 
the  observed  amplification  shows  the  same  frequency  dependence  as  pre¬ 
dicted. 

Another  cw-SACER  experiment  was  conducted  to  determine  the  time 
dependence  of  amplification.  This  experiment  was  similar  to  the  others 
except  a  series  of  tone  bursts  was  sent;  successive  tone  bursts  were 
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TABLE  (4) 

cw-SACER  Data  Comparison 

Mixture:  2-2-16  (H2-C12-SF6) 

Observed  Predicted 

Attenuation  Attenuation  Act 

Frequency  (dB/cm)  (dB/cm)  (dB/cm) 

(kHz)  Before  During  Before  During  Observed  Predicted 


1.0 

.0106 

.0087 

.0096 

.0085 

.0019 

.0011 

2.5 

.0170 

.0156 

.0168 

.0156 

.0014 

.0012 

4.0 

.0240 

.0217 

.0240 

.0223 

.0024 

.0017 

6.5 

.0371 

.0339 

.0386 

.0357 

.0032 

.0029 
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sent  as  closely  spaced  as  possible  without  causing  interference.  A 
1:2:17  ratio  with  a  2.5kHz  tone  burst  was  used  for  this  experiment. 

Table  (5)  and  fig.  (15)  present  the  time  dependence  of  amplification. 

The  point  at  139msec  came  from  a  single  tone  burst  experiment  and  in¬ 
dicates  the  repeatability  of  cw-SACER  experiments.  The  successive  tone 
burst  experiment  can  be  compared  with  an  expression  for  Aat  developed  by 
Toong7  (see  eq.  (55)).  Aat  is  the  difference  between  at  before  and  during 
the  reactions,  where  the  tone  burst  is  written  as  a  decaying  exponential, 
Ao  exp(-att).  The  expression  for  Aat  includes  the  time  derivative  of  the 
temperature;  this  derivative  can  be  calculated  only  for  the  multiple  tone 
burst  experiment  where  a  plot  of  AT  versus  time  can  be  obtained  using 
eq.(67).  Toong's  expression  predicts  Aat  to  be  1.13  sec-1  at  210msec, 
the  observed  value  was  4.64  sec-1.  Although  the  amplification  for  the 
single  tone  burst  experiments  cannot  be  easily  compared  with  the  theory 
developed  by  Toong,  the  frequency  dependence  of  amplification  can  be 
compared  qualitively.  Fig.  (14)  indicates  that  higher  frequencies  are 
amplified  more,  within  the  frequency  range  studied,  as  predicted  by  Gil¬ 
bert's  program.  However,  Ref.  7  states  that  amplification  should  be  in¬ 
dependent  of  frequency  within  this  range.  At  this  time,  no  reason  is 
given  for  the  discrepancy  between  Toong's  predictions,  and  our  data  and 
Gilbert's  computer  program. 

In  addition  to  the  tone  bursts,  unexpected  signals  were  sometimes 
observed.  The  signals  were  of  two  types.  The  first  type  was  similar 
to  the  tone  bursts,  except  had  fewer  cycles.  They  were  approximately 
70  percent  as  large  as  the  tone  bursts,  traveled  with  the  same  velocity. 


TABLE  (5) 


Time  Dependence  of  cw-SACER  Amplification 


Mixture:  1-2-17  (H2-C12-SF6) 
Frequency :  2.5  kHz 

Attenuation  Act 

Time(ms) _ (dB/cm) _ (dB/cm) 


0  100 


200  300  400  500  600 

Time  (ms) 


700 


Figure  15 

Aa  Versus  Time  for  cw-SACER  Experiments 


and  sometimes  decayed  exponentially  with  the  same  attenuation  coefficient. 

At  first  we  suspected  that  the  oscillator  was  being  triggered  by  electri¬ 
cal  noise  associated  with  the  lamps.  However,  the  electrical  noise  hypo¬ 
thesis  was  later  disregarded  after  monitoring  the  oscillator  output  while 
turning  the  lamps  on. 

The  second  type  of  unexpected  signal  appeared  as  a  series  of  pulses 
at  unpatterned  intervals;  each  pulse  was  itself  a  series  of  closed  spaced, 
positive  and  negative  going  spikes.  The  pulses  varied  in  amplitude  from 
about  the  same  size  as  the  tone  bursts  to  over  ten  times  larger;  their 
width  varied  from  about  1  to  10msec.  This  type  of  signal  was  sometimes 
observed  at  low  Cl2  concentrations,  but  occurred  frequently  at  higher  Cl2 
concentrations.  The  unexpected  signals  made  cw-SACER  experiments  at  higher 
Cl2  concentrations  impossible  because  they  completely  masked  the  tone 
bursts.  Repeated  testing  could  not  reproduce  either  type  of  signal  without 
a  reacting  mixture  in  the  tube.  Therefore  we  suspect  that  they  were  both 
true  acoustic  signals.  Neither  type  of  signal  has  been  predicted  or 
experimentally  observed  before  but  could  be  the  result  of  an  instability. 

Attenuation  Measurements 

During  each  attenuation  experiment,  a  tone  burst,  produced  by  the 
sending  transducer,  was  allowed  to  propagate  back  and  forth  in  the  tube 
while  the  amplitude  was  monitored  with  the  receiving  transducer.  The  ampli¬ 
tude  of  the  first  received  tone  burst  was  measured  for  different  transducer 
separations,  giving  the  amplitude  for  various  propagation  distances.  Fig¬ 
ure  (16)  shows  the  results  of  a  typical  attenuation  experiment;  the 
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Figure  16 

Pressure  Amplitude  Versus  Propagation  Distance  for  Attenuation  Experiments 


pressure  amplitude  is  plotted,  on  a  log  scale,  versus  propagation  dis¬ 
tance.  A  straight  line  drawn  through  the  points  was  used  to  calculate 
the  attenuation. 

Attenuation  measurements  were  conducted  for  a  number  of  H2-SF6 
and  CI2-SF6  mixtures,  and  for  different  pressures  and  frequencies. 

Tables  (6-10)  show  the  results  of  these  experiments. 

A  computer  program  developed  by  Shields23  was  used  to  calculate  the 
attenuation  due  to  vibrational  relaxation  by  subtracting  tube  losses  and 
classical  attenuation  from  the  total  attenuation.  The  excess  absorption, 
attributed  to  vibrational  relaxation  was  used  to  determine  vibrational 
relaxation  times  for  each  mixture.  The  log  of  the  relaxation  absorption 
per  wavelength  was  plotted  versus  log(f/p)  and  the  best  straight  line 
with  slope  1  was  used  with  eq.  (60)  to  determine  the  quantity  (px)-1  for 
each  concentration  of  Cl2  and  H2  in  SF6.  Table  (11)  lists  the  vibrational 
relaxation  times;  (px)  is  pressure  times  the  vibrational  relaxation  time 
of  a  given  C12-SF6  or  H2-SF6  mixture,  (px)m  is  pressure  times  the  relaxa¬ 
tion  time  for  SF6  with  only  collisions  with  m  allowed,  see  eq.  (61).  For 
our  experiments,  m  was  either  Cl2  or  H2.  (px)m  should  be  independent 
of  the  concentration  of  species  m,  but,  (px)  varies  with  concentration. 

Additional  measurements  were  taken  to  determine  the  reflectivity  of 
the  transducers  used  in  the  cw  system.  A  tone  burst  was  propagated  back 
and  forth  inside  the  tube  and  monitored  by  the  receiving  transducer.  The 
first  time  the  tone  burst  reached  the  receiver  it  had  propagated  down  the 
tube  once;  the  second  time  the  tone  burst  reached  the  receiver  it  had 
propagated  down  the  tube  three  times  and  had  been  reflected  twice;  etc. 


TABLE  (6) 


Attenuation  Experiment  it  A1 
Gas  Mixture  100%  SF6 


Pressure 

(Torr) 

Frequency 


Attenuation 

(dB/cm) 
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TABLE  (7) 

Attenuation  Experiment  it  AT3 
Gas  Mixture  20%  Cl2/80%  SF6 


Pressure 

(Torr) 

Attenuation 
(dB/ cm) 

Frequency 

(kHz) 

400 

200 

100 

50 

25 

2.5 

.016 

.028 

.024 

.079 

.133 

4 

.035 

.057 

.096 

.157 

.265 

6.5 

.064 

.105 

.205 

.320 

.528 

10 

.139 

.209 

.368 

.694 

1.00 

16 

.319 

.468 

.984 

1.51 

2.28 

25 

.937 

— 

1.86 

2.94 

3.92 

35 

1.51 

3.05  | 

3.38 

7.07 

6.59 

40 

3.45 

5.61 

5.53 

6.65 

6.59 

TABLE  (8) 

Attenuation  Experiment  #  AT6 

Gas  Mixture  50%  Cl2/50%  SF6 

Pressure 

Attenuation 

(Torr) 

(dB/cm) 

Frequency 

.019 

.027 

.05 

.052 

.06 

.107 

.12 

.22 

1"* 

CM 

.47 

.67 

1.26 

63 
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TABLE  (9) 


Attenuation  Experiment  #  AT4 
Gas  Mixture  10%  H2/90%  SF6 


| 

R 


I 

1 

R 


^T»ressure 

(Torr) 

Frequency 

(kHz) 

400 

Attenuation 
(dB/ cm) 

200  100 

50 

25 

2.5 

.011 

.011 

.064 

4 

.024 

.034 

.107 

6.5 

.020 

.046 

.076 

.124 

.197 

10 

.026 

.071 

.124 

.222 

.364 

16 

.058 

.137 

.251 

.440 

1.51 

25 

.142 

.272 

.535 

2.29 

1.56 

35 

.281 

.478 

1.26 

2.09 

40 

.368 

.725 

1.60 

3.48 

3.81 

Mg 


hA  .hi 


.  Jm. ,  fc.i 


TABLE  (10) 


Attenuation  Experiment  #  AT5 
Gas  Mixture  20%  H2/80%  SF6 


Pressure 

(Torr) 

Frequency 


Attenuation 

(dB/cm) 


Inside  the  tube,  the  tone  burst  amplitude  is  given  by  Ao  exp(-ax-nB) 
where  n  is  the  number  of  reflections  the  tone  burst  had  undergone  and  g 
is  the  reflection  coefficient.  Using  a  single  transducer  separation  and 
a  mixture  with  known  attenuation,  the  reflectivity,  exp(-g)  was  calculated 
for  three  pressures  of  pure  SF6j  400,  200,  and  100  torr  (the  frequency 
used  was  4kHz).  For  each  pressure,  the  reflectivity  was  found  to  be  1.0 
within  experimental  error  (less  than  1  percent).  Therefore,  attenuation 
caused  by  reflections  was  not  included  in  absorption  calculations  for 


cw-SACER  experiments. 


CHAPTER  V 


CONCLUSIONS  AND  SUMMARY 

An  experimental  investigation  undertaken  to  observe  sound  amplifica¬ 
tion  in  a  chemically  reacting  Cl2-H2-inert  gas  has  been  described.  A 
pulsed  SACER  was  constructed  for  preliminary  study.  Intense  flash  lamps 
initiated  the  reactions  which  resulted  in  sound  amplification  for  about 
10msec.  A  gain  of  1.8  at  2.5kHz  was  observed  compared  with  a  gain  of  4 
predicted  by  theory.  This  discrepancy  has  not  been  fully  explained  al¬ 
though  it  may  have  been  caused  by  an  error  in  predicting  the  flash  lamp 
intensity. 

Following  the  preliminary  investigation,  a  more  thorough  study  was 
undertaken.  This  study  included  a  broader  range  of  frequencies  and 
initial  conditions-.  During  these  experiments,  low  intensity  fluorescent 
lamps  continuously  dissociated  molecular  chlorine;  the  reactions  progres¬ 
sed  for  a  longer  time,  resulting  in  observed  amplification  for  over 
800m8ec.  The  fluorescent  lamps  were  more  controllable  than  the  flash 
lamps,  so  more  physical  parameters  could  be  varied;  this  made  it  possible 
to  note  trends  in  the  amplification  data. 

For  the  cw-SACER  experiments,  amplification  increased  with  frequency, 
the  values  at  6.5kHz  were  about  2.5  times  larger  than  the  values  at 
1.0kHz.  This  agrees  with  predictions  obtained  with  Gilbert's  numerical 
solution6,  however,  Toong7  states  that  amplification  should  be  indepen¬ 
dent  of  frequency  within  the  range  studied.  The  observed  values  for  Aa 


were  about  20  percent  larger  than  Gilbert's  program  predicted,  but  this 
difference  might  be  explained  by  inaccurate  determination  of  lamp  intensity. 
Toong  has  derived  an  expression  for  Aat7,  the  change  in  the  attenuation 
coefficient  related  to  exponential  decay  with  time,  which  predicts  an 
amplification  only  25  percent  as  large  as  was  observed.  It  is  obvious 
that  the  assumptions  on  which  that  derivation  is  based  are  not  fulfilled 
by  this  experiment. 

A  cw  experiment  was  conducted  which  observed  the  time  dependence 
of  amplification.  'For  the  conditions  chosen  (C12:H2:SF6  ratio  of 
1:2:17  at  l/3atm  and  a  frequency  of  2.5kHz)  amplification  increased 
rapidly  for  the  first  300msec,  peaked  at  about  400msec,  then  started  to 
decrease.  The  temperature  of  the  mixture  followed  a  similar  curve.  As 
the  gas  heats  up  most  of  the  reaction  rates  increase,  so  it  is  not  sur¬ 
prising  that  amplification  increases.  A  point  is  probably  reached  where 
the  supply  of  atomic  H  becomes  depleted  and  reaction  (64),  a  very  exo¬ 
thermic  reaction,  slows  down  to  the  point  where  more  heat  is  lost  through 
the  tube  walls  than  is  generated  by  the  reactions. 

For  the  cw  experiments  in  which  only  one  tone  burst  was  sent  down 
the  tube,  amplification  increased  with  Cl2  concentration.  This  is 
expected  because  increasing  the  Cl2  concentration  also  increased  the 
amount  of  atomic  Cl  available  for  reactions  (64)  through  (65). 

Two  types  of  unexplained  signals  were  observed  during  the  cw  experi¬ 
ments.  Both  appear  to  be  real  since  they  could  not  be  duplicated  without 
a  reacting  mixture  in  the  tube.  The  first  type  was  similar  in  amplitude 
to  the  tone  bursts,  traveled  with  the  same  speed,  and  sometimes  decayed 
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exponentially  with  the  same  attenuation  coefficient.  The  second  type 
was  much  larger  than  the  tone  bursts  and  did  not  propagate  back  and 
forth  the  tube.  The  second  type  was  observed  frequently  when  the  CI2 
concentration  was  large  or  when  the  temperature  was  high. 

It  is  clear  that  future  studies  of  amplification  in  reacting  systems 
is  merited.  Both  of  our  studies  were  hampered  by  the  question  of  light 
intensity  within  the  tube.  Perhaps  a  fiber  optics  probe  could  be  used 
to  determine  light  intensity  within  the  tube.  Some  measurements  have 
been  made  in  our  lab  with  a  fiber  optic  probe  positioned  outside  the 
tube,  but  with  little  success.  Two  types  of  unexplained  signals  were 
observed  during  the  cw  experiments.  Since  neither  signal  was  theoreti¬ 
cally  predicted,  a  set  of  experiments  should  be  conducted  to  determine 
their  origin.  A  type  of  chemical  instability  may  be  responsible.  A 
final  suggestion  for  further  study  is  to  develop  a  system  with  flowing 
gases  so  that  the  reactants  are  continuously  renewed.  Such  a  system 
would  allow  signal  averaging  and  thus  lessen  experimental  error. 
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Program  RICABS.F4 


This  program  was  used  with  Gilbert's  computer  program6  to  predict 
sound  amplitude  for  cw-SACER  experiments.  A  similar  program  was  used  for 
the  pulsed  SACER  predictions. 
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II.PLICIT  pf;aL*3  CA-H,0-Z) 

DIMENSION  XYl  J),Pkr(3,3),V.Z(3),WiS(3) 

DIMENSION  ARX(  3 ) » AIX ( 3 ) »  BFOR ( 3 ) , BFOI ( 3 ) 

DIMENSION  PBP(3),RBIC3),BF  lf.(3),BFlI(3),UT(b) 

DIMENSION  T Y ( 30 ) 

CC 

DIMENSION  Y(79),m(f9),Y12(79),Yi3t79),  Y22C»9) 

CCC 

SURTKIX ) sDSuRT (X) 

C0SMX)*DC03CX) 

SINF(X)«DSIN(X) 

EXPF(X)*DEXP(X) 

CC 

N»79 

DO  J*l,N 

READ (10# 4 JO)  Y(J), Y1 1  ( J) ,  Y12 ( J> , Y1 3 ( J) , Y22 CU ) 

400  CONTINUE 
410  FORMAI(SF) 

CCC 

1  FQRMAI ( 5D/4P) 

2  FORMAT (  2D ) 

100  FORMAT ( 6H  El V*D1 3 .6# 1X6H  E1A-D13.6) 

102  F0RMAK4H  T»D12.5,IX4H  W«D12.5,IX5H  C8=Di2.5, 1X5h  CV«bl3.6) 

103  FORMAT ( 7H  VIS  *012.5, 7H  P (MM)*C 1 2. 5 , 7H  R (CM ) *D12 .5, 4H  FM=D12.5) 

104  FURMAT(7H  TA+CA=DI2.5,4H  RA*D12.5, 1X6HWLCM)=D12.5, 5H  FSVSD12.5) 

105  FORMATISH  ACS*Dl2.1,oH  ACT*Dl 2 . 4, 3X, ' SRTMG= ' , 01 2 .5 , // ) 

106  FORMAT ( 5H  RAP=DI2.5,6H  FSVP=01 2 . 5, 1X6HALFAP»P12 . 5, 6H  VELP=D12.5) 

«07  F0RKAT(4H  RW*012.5,5H  RP*D12.5,7H  GAMAK*D1 2 . 5 , 7H  GAMAI=D13.6) 

102  FORMAT ( IX, 'F«' ,i 5. 1 , 2X, ' MS* ' ,09. 3, 2X, * V« • , DIO. 4, 

*2X, •V(%)*',F4.0,2X, 'V/W=' , D10.3, •  /',F6.0,/) 

Wl*146. 

W2*7 1 , 0 
W3*2.o 
R«1.23 
C  R*1 . 27 

ACS*1 . 

ACT*!. 

CC 

SM1*5 . 5 1 
SM2»4.*0 

SM3*2.9 

EOK3*37 . 

EOK1*200.9 
EOK2»257.o 
SM12*f  5M1  tSM2)  /2. 

Sh23*(SM2*SM3) /2. 

SMl3*(SMl+SM3)/2. 

EOK12*(FOK1*EuK2)**0.5 

LOK23«(EOK2*EOK3)**0,5 

E0K13«(E0Kl«e0K3)**0.5 

TABS*0. 

TDL«0. 

READ (7,665)  FREU 
F«FREQ/<2.*3.I4I6)/I000. 

WRITE (8,880)  F 

880  FORMAT! 2X, 'fREQa',F10.2,2X, 'KHZ.',/) 

WRITE! 8,890) 

890  FORMAT! 3X, 'TIME' ,6X' TEMP' ,4X, 'P(ATM) ' , 2X ' AMP ' , 4X» ' NEW AMP ' , 

12X' TOTAL APS ' , 2X ' DIST ' , 5X ' DB/M ' , /) 

66$  F0RMAT(9D) 


CV1 *3 • 0 
CV2«2.S 
CV3«2.5 
C81»t.5 
C82«1.5 
C83«1.5 

READ(7, *65)  TIME, TEMP,  PATM, AMP, FM1 ,  f  M2 » FM3 
AMPNEW»AMP 

WRITEC  8,  800)  TIME, TEMP, PATM,  AMP,  At'PNEW , TABS 
800  FORMAT (1X,F3.5,F9a2, 3F8.3,4f9.5) 

>00  CONTINUE 

OTIME»TTME 

REA0(7,665,£ND»99)  TIME, TEMP, PATM, AMP, FM 1 ,  FM2, FM3 

p«760.*PATM 

T*TEM° 

TS1«T/E0K1 
TS2»T/EOK2 
TS3«T/EOK3 
TS23*T/EOK23 
TS13*T/E()K  1  3 
T512»T/F0K12 

CALL  SSEVALCN, TS12,Y,  YU. SEVAL) 

0M11«SEVAL 

CALL  SSEVALCN, T312,Y,Y12»3EVAL) 

0MI2»AEVAL 

CALL  SSEVALCN, TS12.Y.Y13. SEVAL) 

0M13«SEVAL 

CALL  SSEVALCN, TS12,Y,Y22»SEYAL) 

OM22-SEVAL 

CALL  SSEVALCN,TS1,Y,Y22, SEVAL) 

OM221«SEVAL 

CALL  SSEVALCN, TS2,Y,Y22, SEVAL) 

OM222»SEVAL 

CALL  PDO»1,W2,TEMP,SMI2,OMU,PRU(X,2),PRD(2, 1)3 
TYPE  <»27,FRDU,2) 

)7  FORMAT  C  2X, 1  PHD (1,2)*' ,  D 1 2 . 5 ) 

VVIS(l)«20o.93*10**-7.*(Wi»T)**0,5/(5Ml*SMl*0M221) 
VVIS(2)»206.93*l0**-7,*(W2*T)**0.5/CSM2*5r.2«OM222j 
TYPE  937, Wlj(2) 

>7  FORMAT C2X, ' VVI5C2)»'»D12.5) 

•  TYPE  *02,TS1,TS2,T512,(JM11,OM221,OM222 

402  FORMAT (2X,6F9,3) 

FM2R1.-EM1-FM3 

XA»(FM1/W1) +CFM2/W2)*CFM3/W3) 

XJ«FM1/W1/XX 
X2*FM2/W2/XX 
X3>FM3/W3/XX 
M*X1*i»1*X2*W27X3#W3 
TYPE  ^47,X1,X2,X3,FMJ,FK2,FM3 
17  FORMAT(?X,6(2X,C12.5)3 

CV«Xl*CVl7X2«CV27X3*CV3 
C8»X1*C81+X2*C82+X3*C83 
CALL  SSEVAL  CM , TS23, Y, Y1 1 , 5EVAL) 

0M1 IsAEVAL 

CALL  SSEVALCN, TS23 , Y, Yl 2, 5EYAL ) 

0M12»SEVAL 

CALL  SSEVALCN, TS23,Y,Y13, SEVAL) 

0M13-5EVAL 

CALL  SSEVALCN, TS23,Y,Y22,SEVAL) 


v'-S' v. 


OM22=*EVAL 

CALL  SSEVALIN,TS2, Y,Y22,SLVAL) 

0M22 1 *SEVAL 

CALL  SSEVALCN,To3,Y,Y22,SEVAL) 

0M222»SEVAL 

CALL  PD(W2,w3,TEMP,SM23,omi,PRD(2,3),PRD(3,2)) 
WlS<3)*26o.93*10**-7.«<W3*T)**0.5/(SM3*SM3*QM222) 

CALL  SSFVAL(N,TS13,Y,Y11»SKVAL) 

OMl laSEVAL 

CALL  SSEVAL(K,TS13,Y«Y12,S£VAL) 

OMliaSEVAL 

CALL  SSEVALCN,TS13,Y, Y1 3,5EVAL) 

OM13*SEVAL 

CALL  SSEVAL(N,TS13,Y,Y22,SEVAL) 

OM22*SEVAL 

CALL  5SEVAL (N,  TSl,Y,Y22,Sfc.VAL) 

OM221»5EVAL 

CALL  $SEVAL(N,TS3,Y, Y22,SLVAL) 

QM222»SEVAL 

CALL  PU(wi,Wj,TEh.P,SKl3,OMll,PRD(l, 3) ,PRDC3, 1)3 
CACULATION  UF  VISCOSITY  FUR  3  GASES  USING 
DIFFUSION  CONSTANTS  FOR  TWO  GASES 
XY(1)*X1 
XY(2)«X2 
XY ( 3) *X  3 
WZ(1J««1 
WZ(2)«W2 
KZ<3)«*3 
VIS*0. 

DO  990  I«l,  3 
PPD»0 

DO  980  K»l,3 
ir(K.EO.I)  GO  TO  980 

1PD»PP0  ♦  C  t .  385*XY ( I)*XYUO  *8,314403/ (PHD ( I, K)*WZ( I ) ) ) 
FORM AT  f  4D 1 2 . 5 ) 

TXPE  933, PPD 

TYPE  9l7,XY(I),XY(K),PPDCI,iO,WZCI) 

CONTINUE 

IF(XY(I).LT.1D-5)G0  TU  990 
VI5»VIS  ♦  (XY(I)**2./CXY(I)**2./WISCI)*PPD)) 

TYPE  913, VIS, VVIS(I), PPD 

FORMAT(2X, ’VIS»'#Dt2.5, 'VVIS»',D12.5, ’PPD-' ,D12.5) 
CONTINUE 

TYPE  933, PPD 
FORMAT<2X,'PPD»',D12.5) 

TYPE  939, VIS 

FORMAT (2X, 'VIS*1 ,012.5) 

CONTINUE 

V»VIS 

RTM«8.3144D7«T/W 
ZROT«83.3/EXP(16.7/T»*. 333333) 

FM»(2 ,5*1 . )*1 .25*101 4000, /V/ZROT/ (6, 2832*2 ,5) /I . 
SRTMG»SQRTF(RTM*< 1 ,D0+1 .DO/CV) ) 

WRITE(3,102)T,W,C8,CV 
WRITE?  3, 103)V,P,R,FM 
WRITE f  3 » 1 05 )  ACS, ACT,SRTMG 
AF«2.D0*3. 1416D0*1000,O0*F 

AT«F*760,D0*SURTF<CV*(CV*1. DO) /(C8*(C8*1. DO) ))/(FM*P> 
TYPE  973, FM 
FORMAT (2X, 'FM«', 012,5) 
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WX,*SRTMG/(F*1000.)/100. 

PDN*l.3332C3*P 
RTM»8 • 3 1 4407*1 /W 
YOD«V*RTM/PDN 

Cl«Xl*t.S*X2*1.0*X3*1.0 

C2*7.f92 

AT1»AT 

PTSF6*7.25E-7 
PTCL2«4.37E-7 
t'TH2»l «7d£*8 

PTMIXal  ./(Xl/rTi>F6+X2/PTCL2*X3/PTH2) 
AT2»7t>0./F*pTMIX 
C  ATl«1.0E-9 

C  AT2»1.0F-9 

C  TYPE  4(>7 »  ATI  i  AT2 

967  F0R«AT(2X, 'ATI*' ,D12.5)2X, ' AT2» • ,01 2.6) 

AF2*AF*AF 
PAt»t .♦AT1*AT1 
PA2»t.+AF2#AT2«AT2 
Cl»CV-C» 

C2»X1*C2 

A2*CV*CV+C8#C8*AT*AT 
CVP«CR*C1/PA1+C2/PA2 
CVI«“(C1*AF*AT1/PA1+c2*AF*AT2/PA2) 
CVAV2«CVR*CVR+CVI*CV1 
UFN«VUD*<1.00*2.25C0«CVf'./CVAV2) 
UINa-2.25O0*CVI/CVAV2*VQD 
RSHR«1.D0*CVR/CYAV2 
P5HI«-CVI/CVAV2 
8R»-RTM*R5Hfl*AF#UlN 
BX»-RTv*rshi-AF*(URN*1.3333300*VOD1 
AP.*3.33J3DO*VuD*URN+RTM*U1N/AF 
AX«l.3333f>0*VOU*UIN-RTM*URN/AF 
AB«2»00*t  AR*AR+AI*AI) 
ABR«»(8R»AR't6I*AI)>'A8 
Abl»-(8I*AR-BR*AI}/AB 
Ar,AR«AR*AP  +  Al*AX 

FF.«ABR*aHR-A81»A8X*AF*AF*aP/ARAR 
ri»2.00*A8P*ABI-At *AF»AI/APAP 
FRI»SORTF  CFR*FR+FI*FI) 

ircru  it)i2(i2 

11  GR*SQRTF(*(FPX+FR) '2,U0) 

Go  TO  13 

12  GR«5QRTF((FRI*FK)/2.C0) 

13  GI«8QRTFt(PRI-FM/2.0D0) 

RL2»ABR*GR 

FI2»A»I*GI 

P.U«ABR-GP 

FI1»ABI-G1 

AT*  •P/'7t>O.DO 

RF»5RTMG/AF 

RFS«RF*RF 

GISR«R1)1*HFS 

G1SI»FU*PFS 

G25R«R1<2*RFS 

G2SI»FI2*RFS 

RV»R/RF 

C*  TYPE  501 (GISR)GISX 

C*  TYPE  SOI )G2&R)G2SX 

SOI  FORMAT ( SO /4D) 
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RViaSRTMGaSRTMG/  (  VuD*AF) 

SglGRaHRtJ/ VOD 

SQIGXal'IK/VOD 

GlSAVaG1SR*GlSR+ClSI*GlSI 

BETlRaGlSWGlSAV-tiGICR/RW 

BETlI«G>SR/GlSAV-5QIGI/kW 

G2SAV*G2SP*G25R+G2SI*G2SI 

8LT2P»G2SI/G2SAV-SQXGR/RW 

BLT2l*G?SP/G2SAV-5QXGI/RW 

FKPAaFMTM 

FKaFKPA'FM 

G2ERO»1.0n+(t.00/CV) 

GINF*1.DO*(1.1,0/C8) 

CGaGZERG*(l.O0+FN*FN)/(l.P0+(GZER(j*FN*FN/GIHF)) 

TC*(CG-1 .DO) /CG 

Tb»C9.DO*CG-5.D0)/4.D0 

TKal  3 . 3-)3900«(  1  .DO/SOF.TFtCG)  ♦TC*SGRTF  (Tb)  )*5QRTF ( V )  /R 
TAaTK*5(JRTF(F/P) 

FaVaSGRTp ( R .31 44D3«CC*T/W ) 

C  TXPE  940, FSV 

940  FORMAT ( 5F ) 

DELV«TK*FSV*FSV*.0018323DO/SQRTFCF*P) 

CA«( 1 .3  33  33300+TC*TB)«V*SQRTF (H/T) 
CAal4.106D0*CA*f«F/{P*S0PTF(CG*CG*CG)) 
RAa3l4l.62DO*{CV-C8j/iORTK(CV*l.DO)#CC8»l.DO)*CV*C8) 
VEL«(F5V-OELV)«100.D0 
C  TFPE  «*40, VEL 

TVaVEL 

RA»RA*fM*F / <VEL*t 1 ,DOtFH*FH) ) 

ALPHA»RAMTA+CA)/8. 68t,00 
TCA«(TA*CA)*100. 

ALFA* ALPHA 

PCRstALPHA 

PCI»AF/VEL 

PC2R»PCP*PCR-PCX*PCI 

PC2I«2.P0*PCP*PCI 

C*  WRIIEO.104)  TCA,RA,HL,F5V 

7  ZH«PC2T-»F/VgO 

ARARaSQBTF (PC2R*PC2R+21I*2lI) 

ARlRa-R*SORir C (PC2P+AKAR)*.5DO) 
ARH«R*5Q»TF({-PC2PtARAR)*.5D0) 

Z3R»PC2B-PL2 

Z3I«PC2I-FI2 

AP.AR»SQBTF(Z3R*23R+23I*23X) 

AR3Ra*R*S0RTF( .SDQ*( ARAR+Z3R) ) 

AP.3IaR*8QBTF(  ,5D0*C  ARAP-Z3P)  ) 

22R*PC2«-PL1 

22l»PC2I-FIl 

AHARa5ttHTF(Z2R*Z2K+Z2X*Z2I) 

AR2R>R*8CiRTF(  .  5D0*  ( ARAR+Z2R)  ) 

AR2l«R*«QRTF( ,5D0*( ARAR-Z2R) ) 

ARX(  1 )  ■  AMR 
ARX(2)«  AP2R 
AKXO)*AR3R 
AXX(l)a  API I 
AIX(2)a  AP2X 
AIX(3)»AR3I 
DO  30  J»t«  3 
ARa  AKX(J) 

Ala  AIXlJ) 


AA*5rtf'TFCAP*AF  +  Al*An 
If CAA-6.319, 19,25 

19  CUSX«AR/AA 
SZNXaAI/AA 

CUS2X»l.DO-2.U0*SINX*SINX 

SIN2Xs2.D0«C05X*S1NX 

B1I«AA*SINX/2.D0 

BlR«AA*COSX/2.D0 

FlMsAA/J.00 

DU  20  I»l,30 

X*I 

COSPX»CDSX 

COSX»COSX*COS2X-SINX*SIN2X 

SINX«»INX«CQS2X+CUSrX*SIN2X 

FlM«-n«*AA*AA/(4.P0*X*lX+l.D0)) 

81 lRaCOSX*FlM 

61RaBlR*BUR 

BUIaSINX*FlM 

B1I«B1I*BUI 

If (FI. 00000001)21, 21,20 

20  CONTINUE 

21  CUSXaAR/AA 
SINX«AI/A» 

CUS2X«l.DO-2.D0*SINX*SINX 

SIN2X»2 ,DO*COSX*SINX 

CC5X»C082X 

SINX»!»IN2X 

FOM«-AA*AA/4.DO 

BOR»l ,DO-AA»AA#CO6X/4.D0 

B0I»-AA*AA*5INX/4.D0 

DU  22  I»2,30 

X*I 

FOM*“FOm*AA*AA/ l4,PO*X*X3 
CUSXP»COSX 

CUSX«C0SX*C0S2X-SINX*5IN2X 

SINX«3INX*C0S2X4CUSXP*SIN2X 

BOlR»FO*«COSX 

B01I«FO***SIMX 

BOP»BOR+BO  1R 
BOI«BOI*BOU 

IF (F0M*F0M«, 0000000 1)23, 22,22 

22  CUNTINUF 

23  GO  TO  20 

25  AR«-AR 

Ala-Al 

SRZRaSQPTFC  ,5D0*(  AA+AR) ) 
3R2I*-50RTF ( . 5D0* 1 AA-AR 3 3 

18  AA2«A«*AA 
AA3*AA2*AA 
AA4«A*2*AA2 
AA5«AA4*AA 
AA6«AA4*AA2 

ClR«.797fl846D0*SRZR/AA 

C1X«-.7978846D0*SRZI/AA 

AR«R«AR*AR-AI*AI 

AR8R>ARSR/AA2 

Alt8I»2.O0*AR*AI 

AR5I-ARSI/AA2 

ARCR»AR*ARSR-AI*ARSI 

ARCRaARCR/AA 
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AkCI*AI*APSR+AR«AkSI 

AKCJ»APCI/AA 

AFFRsAF'SR*ARSP«ARiiI*ARSI 

ARFk*ARFR/AA 

ARFI»2.00*APSR*ARSI 

ARFl»APEl/AA 

AnPR*ARFR*AR-APFI*AI 

ARPR«*RPR/AA 

ARPX*ARpP*AI+ARFI*AR 

ARP1*APPI/AA 

arhr»arcr*arcr-arci*akci 

ARHftaARHP./ AA 
AKHI*2.D0»AFCR«APCI 
ARH1»APHI/AA 

CERlsl  .I>0*.117ia8uO*ARSP./AA2-.14420DO*AP.FP./AA4 

CBRl«ChPlA.67659279D0*ARHR/AAo 

CBI1«-. t 1 7188D0*ARSI/AA2+ . 1 4420l)0*AKFI /AA4 

CBIl«CBl 1-. 67*5927900* ARHI/AAb 

SBRi«.375D0*AR/AA2-.10254D0*AKCR/AA3*.2775ai;0*ARPK/AA5 
SBI1«-.375D0*AI/AA2*.10254D0*ARCX/AA3-.2775BD0*ARPI/AA5 
Al»COSF( AP-2. 35b 195&0} 

A2»SINFfAP-2.356l95D0) 

IF ( AI .LT . 40 . )  GO  TO  201 

EPal.OQO 

EMsO.ODO 

201  IFUI.GT.-40.}  GO  TU  202 
EP*O.ODO 

EUal.ODO 
GO  TO  203 

202  EP*EXPF ( AI ) 

ENaEXPF ( -AI ) 

203  CONTINUE 
CFR1»A1«(FP*EN)/2.P0 
CFI1»-A2*(EP-EN)/2.D0 
SFR1»A2*( EP+Eh ) /2.D0 
5FI1«A1*(EP-EN)/2.P0 

B»C1F*CHP1*CFM-C1P*CBI1*CFU-C1I*CBP1*CFI1-C1I*CBI1#CFR1 
BlP«B*CiP.*SBRl*£FRltClR*SbIl*5m*ClI*SBRl*tFlUClI*SfcIl«5FR1 
BI»C1I*CBP 1*CFR 1 tC 1P.#GBI1*CFR ltClk*CBkl*Cl I 1"C 1  I*CBI 1*CFI l 
81X«B1-C1  t*SBRl#$FRl-ClK*5BIl*SFPi-ClR*SBFU*SFIl*CiI*SBIl*SUl 
CBR0»l.O0-.0703125D0*APSR/AA2t.U2152DO*AhFR/AA4 
CBR0«CBP0-.57250U>0*ARHR/AA6 
CBI0«t.o7O3125D0*ARSI/AA2-.112152L0*ARri/AA4 
CEl0«CBI0*.57250U/0*ARHl/AAb 

SBP0«».125P0*AP/AA2+.073242D0*ARCR/AA3-.227J08D0«APPP/AA5 
SBI0«, 12500a AI/AA2-. 073242D0*ARCI /AA3+. 227 108D0*ARPI/AA5 
A5«C0Sr(AP-. 785399200) 

ABaSINFCAR-. 785398200} 

CFF0«A5*(EP*Eh)/2»P0 
CFlOa-AF#(EP-EN) /2.D0 
3FR0«A6»(EP*Fh)/2,D0 
Sri0«A5*(EP-EN)/2.D0 

B0R»ClR*C8R0*CrR0-ClR*CBI0*CFI0-ClI*CBR0*CFI0-ClI*CBI0*CFR0 

BOR*BOR"ClR*3BRO*SFROTClR*SBXO*sriO*ClI*S8RO*SFIO*Cl I*SBIO*SFFO 

BOI«-Ci:*CBXO*CFIOaClR*CBRO*CFIOatlR«CBIO*CFROaCtX*CBROaCFPO 

BOI*BO:*ClI*SBIO*5riO-ClR*BBPO*5FIO-ClP*SBIO»SFRO-ClIaSBPO*6FRO 

B1P«-81P 

BII»-B1I 

26  BF1B(J)«B1P 
BFll(J)»BH 


*m*  **  *  *  n  -'1 


BF0R(J)»BOP 

BF01(J)*B0I 

C**  TTPE  '>99,J,blR,bll,bOR»BOX 

Dl*flFOB(J)*ftFOP(J)+iFOrCJ)*UFOICJ) 
RBR(J)*<BF1R{J)*BF0R  (J)*BFH(J)*BFOI(J))/bl 

30  RBI(J)*(BF1I(iJ)*BF0R{iJ)-BF1R(J)*BF0ICJ))/DI 
DEL*SaRTF<  3. 1410DO*GZERO*.5DO) /(RP*RW) 
DELS«DEt.*(2.D0-ACS)/ACS 
DEbtauEL*4.5D0*( 2. DO-ACT) /ACT 
DELTD»( SQ I GR+ 3 . 5D0 ) *(SQIGR*3 , 5DQ ) tSQICI*SQIGl 
DELTR«DELT*( SUIGR*  CSQIGR+3 . 5D0) +SiiIGI*SuXGI ) /DELTD 
DELTI»DF.LT*(SUIGI*C bQIGR+3 . 5D0 ) -SQIGR*SUIGI ) /DELTD 
01 1R»1 «PO-DELS*l AR1K*RBR( 1  )-ARH*RBX  ( 1 ) ) 
DUI«-DELS*CAR1I*RPR(1)+AR1R*RBI(1)) 
Sfcl2P»A»3fi*R8R(3)-AR3X*RBIf3) 
SD12I«A»3p*RBI(3)+AR3I*RBR(3) 
DX2R«l.D0-(DELTR«5D12R-DELTl*5Di2I) 

012Ia-D£LTR*S01 2I-PELTI*SD 1 2R 

AP,2AV*AP2R*AR2R+Ak21*AR2I 

SC 13R»( PET1R*KBK( 2)-BETlI*PtI C2 ) )*RF/AR2AV 

SDX3I»C^ET1R*kBI<2)+BET1I*RbRC2))*RP/AR2AV 

D13R»SD1  3R*AR2Pt5D 1 3I*AR2I 

D13I»-SD1 3P*AR2I*5D13I*AR2R 

SblRaDl  1R*M2R-D1 1X*D12I 

SD1I»D11R*D12I+D1  1I*D12P. 

D1R«S01P*D13R-SX-1X*D13I 
DlI»SDtR*D13I+SDU*D13R 
D21Ral .DO-LELS*( AR3R*RBR<  3) -AR3I*KBI C  3) ) 

021  Xa-DFL3*( AR3R*ftfi  1(3)  *AR  3l*RBR  {  3  )  ) 

SD22R»AP2P*RBRC  2) “AR2I*RBI ( 2) 

SC22I«AB2P*RBI(2)+AR2X*RBR(2) 

D22R«l.DO-(DELTR*SD22R-0EETX*SD22I) 

D23Ia-DELTR*SD22I-PELTI*SD22R 

AR1AV»AP1P*AR1R+AR1I*AR1I 

SC23R«( »ET2R*RBR{  X )  -B*;t2I*PBI  1 1 )  ) *RP/AR  1 AV 

5D23I»(BET2R*f<BX ( l)-BEl2I*PbR(l ))*PP/AR1AV 

D23R»SD23P*AR1R+SLi23I*ARH 

D231«-SD23R*AR1I+5D23I*ARIR 

SD2R»022R*D23R-D22I*P23I 

Sb2IaD22R*D23l+D22I*D23R 

D2R«S02P*021R-S02i*U21I 

02I>S02I«021RtSl>2R*021X 

D31R*  t . DO-DELS* ( AR?R*RBR( 2)-AR2I«RBI(2)) 
C3lI»-DELS«(AP2R*KBI( 2 ) +AK2X*RBR( 2) ) 

Sb 3 2R»AR 3 R*RBR{ 3 ) -AR3I*RRI ( 3) 

Sb32l«AP3i*RBR(3) tAR3R*RBI( 3) 

D32R«l.DO-(OELTR*5D32R-DFLTI*SD32I) 

D32I«-DELTI*SD32R-DELTR*St 321 

5b33R«(RETiR*RBRa)-BETlI*PBICl))*RP/APXAV 

S033I«(HET1R*RBI(1)+BLT1I*PBR(1))*PP/AR1AV 

D33R«AR! P*SD33P+AR1I*SD33I 

D33I»AR1R»SD33I-AR1I*5D33R 

Sb3R>U33R*D32K-D 331*032 I 

8C3X»033R*032I+D33I*C32R 

D3R«S03R*03XR-Sb3X*C31I 

D3laS03P*031I*SP3X*D3XP 

DK«01R*O2R-D3h 

DX«01I*02I-D3I 

E1MD1R 

EXXaOXI 


K21h  =  wl  IP 
£21 1  =  l  1  1  I 
Ei2t-su2>l 
£221*^2  21 

St23Ra(AP3P*Pi;V!(3)-AR3l«PLlC3)  5 /FT 
Sfc23ia(AR3I*At  PC  3)*AP3P*PL  1(3)  )/Pt 
fc.*3fcsSE/3P*r.CT21  -uF23I*h£.T?l 
1 2  3 1  *  a  E  2  3  B  •  r.  F '»  21+^f  23i*btT2F 
J>L2t,*c23t'*£22F”E7JT«F22I 
5t21*£2  3P»E221  +  t.  23I*t22h 
E2PsE21P«S£2P.-E21i*0E2I 
F2I»E2lR*SF21*F.2U*iE2R 
El  sElF*(;iSp.-ElI*01Sl*L2I. 

E1«E  l  A*Ma  T  ♦£  1  I*GlS»\*t.21 

G6ps(f:H»rp*Li*t'i)/(rp*ni-.+Lt»pn 

05IsCEI»DP-t:H»L'I)/(yR*DP.+l'T»DI) 

PC2F RsGSP/RFS 
PC2P1«CM/RFS 

AFAMSQ-TF  (  PC  ipp,*pr  2PK+PC  ^Pl*FC2P  X  ) 

P(.PK«a>JPTF(  ,5LiO*(Pr2PF.tAHAP)  ) 

PCPl«au«TK .5uO*(-PC2PR+»PAH) ) 

VELpaAF/PCFI 
ALFAPapCPR 
EIVa(VELP-VKL)  /VEX, 

EIA«(ALFAP-ALPHA)/ALrHA 

VEL*VELP 

ALPHASALFAP 

PC2R«PC2PP 

PC2I»?C2PT 

IF  (EIV*F'IV-, 00000001)41,40, 40 
Gu  TO  7 

IF (El A*FI»-. 000025) 43, 42, 42 
G(j  TO  7 

PCVP«(VELP-TV)/(DELV*100.) 

GAMAIaSBTMG/VELp 

GAMAR«ALF»P*SRTHG/AF 

RAP*(  AI«FAP“(TA*CA) ^8.686)*(2./F*1000.) * ( VELP+CELV* 1 00 . ) 
ALFAP«ALFAP*8. 666*100. 

ROT*ALF»F-TCA 

Al>FPW«DFA*F*  1000.  *760. /P 

ADF«AOFPW/WL 

ACLDr»ALFAP+ADr 

VAPM«OBPM-ACLBF 

V0AP«VAPH*100,/ACLDF 

VAPa«VAPM*hL 

Dfc»(TIMF-UTIME)*VELP/100. 

TDL«TOt,+Dt, 

DAB5*ALFAP*0L 
TABSaT ABS+ DAPS 
AMPHEW«AKP*10.**C -TABS/20.) 

WRITE!  8, 800)  TIME, TEMP, PATH, AMP, AMPNEW, TABS, TDL,ALFAP 

WRITE! 3,10b)  RAP, F3VP, ALFAP, VELP 
WRITEC 3,107,  PW,RP,GAMAR,GAMAI 
WRITE!  3,  308)  F.l/BPM,  ACLDF,  ADF 

FORMAT! ’  Fa',FS.l,3X, '  MS*' , D9 .3, 3X, 'CL+DIF*',D9.3,3X, 'DIF»' ,t>9.3) 
WPIIE(  3,  108)  F*DBPM»  VELP, VOAP, VAPW ,  FPA 
GO  TO  «»0 
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c*** 
c #** 

SUBROUTINE  VISDF(Xl,X2»R»Wl,W2*T*CV,VIS,DFA,Ot111,QM12»QM13, 
1CIF.22,0M221,CIM222,&M1,SK2,3M12,T31,TS2,TS12) 

IMPLICIT  REAL*8  <A-H,0-Z) 

W12«2.*W1«W2/(W1»W2) 

HM«X1*v<UX2*H2 
GM»1.»1./CV 
Rft*8 .31  44D7 
C«(GM*RP*T/WM)**.5 
PATM«P/760. 

C*** 

C*** 

C*  TIPS  12, 1*3 12 

12  FUHMAr<2X, 'TS12s' ,F7.3) 

20  FORMAT (2X»4F9.5) 

AS  1 2«OM22/OM 1 1 
BS12»(5.*OM12-4.*UM13)/CjM11 
CS12«>JM12/0M11 
SOT  1  2aSM  1  2*SM  1 2*Otl22 
C*  TYPE  1 3,TS1 , T32 

13  F0RMAI(2X. 'TS1.TS2  5  * , 2fS . 3, /, 2X»  ’  0MACA<2,2}1  ,  OMACA ( 2, 2) 2 1 ) 

C*  TYPE  20 , 17M221 ,  OM222 

SOTl*SM1«SMl*UM22l 
SOT2«SM2*SM2*<JM222 
Vl«.0O0ii26693*(V,l*T)**.5/S0Tl 
V2*.0O0026693*CW2*T)*#.5/5OT2 
C  -  CALCULATES  VISCOSITY  OF  THE  MIXTURE  - 

V12«.000026S93*(W12*T)**.5/50T12 
XV»X1**2./V1+2*XI*X2/V12*X2**2./V2 
YV«X1**2.  'V1»W1/W2+X1*X2/V12*(W1+W2) /W12*V1 2**2 . / ( V1*V2 ) 
YV«YV*X2**2./Y2*W2/Wl 
YV».6*AS12*YV 

ZV«Xl**2.*Wi/w2+2,*Xl*x2*(CWltW2>/C2.*W12)*lV12/Vi*Y12/V2)-l.) 

2Y«ZV*X2**2.«W2/W1 

ZV».6*A512*ZV 

VIS«(1*ZV)/1XV+XV) 

C*  TYPE  14, T, VI, V2, VIS 

14  FORMAT! 3X, ’ T!K) 1 » 5Xi • YI5 U ) ' » 5X, ' VIS ( 2) ' , 5X, ' YIS ' 

•#/»  2X#F7«2#3Fll .7/) 

C  -  CALCULATES  DIFFUSION  TERM  - 

PDl2».002628*(T**3./W12)**.5/lSM12*SM12*OI-,lt) 

C 

RL1«,00019891*(T/W1)**.5/SOT1 
RL2«. 0001989 1*!TKW2)**.5/S0T2 
RL12».00019891*(T/W12)««.5/SOTI2 
C 

UA«4,*A31 2/15, 

Ut»»(l2.  *6512/5.  tl.)/l2. 
gC«(Hl-«2)*»2./!Wi*H2) 

C 

Ul«UA-UB*Wi/W2*UC/2. 

U2«UA-UB**2/WI  t  UC/2. 
UY»UA«(Ml*W2)/'(,5*ta2)*RLl2**2,XCRLl*RL2) 
UY«UY-UR-5./!32.*ASl2)*!2.4*BS12-5.)*UC 
C 

XL»X1»X1 /RL1  tZ.PXlPXZ/RLiZ  ♦X2*X2/RL2 
YL»X1*X1*U1/RL1  *2.*X1*X2*UY/RL12  *X2*X2*U2/RL2 
S1«!W1 +*2) /!2.*W2)*RL12/RL1  -(V.2-W1  J /( 2  ,«W  U /UA  -l. 


84 


S2*(N1+rf2)*RL12/(2.*Wl«f L2)  -( WI-W2 )  /  C  2  ,*W2 ) /UA  -1. 

C 

RKT*Xl*X2/(b.*RU2)*(Sl*Xl-S2*X2)/(XL+YL)*C6.*CS12-5.) 

C  TYPE  17,C,GM,1.KT 

17  FURMAT(2X, 'C,GM,RKT*  ’#3011.3/) 

DF*X1*X2*GM*GM*PD12/C**2.*<(W1-W2)/WM  ♦  ( 1 . -1  ,/GM)*PKT/ (XI *X2 ) ) **2 . 
Of A*8.68b*2.*3.1416**2./GK*DF 

AVIS*N.686*2.*3.141b**2»/GM*(4./3,)*VIS*lQ00./ (PATM*  t013000.) 
ATHM*AVIS*.75*(l.+2.25/(.V)*(GM-l.)/GM 
ADFPK*DFA* 1000. *760. /P 
C  TYPE  13.  AVIS, ATHH.ADFPK 

18  FORMATC2X, 'AVIS,  ATHM  t  ADFPK(Db/WL/KhZ)  3011.3///) 

RETURN 

END 

C*** 

SUBROUTINE  SSEVAL(N,U,Y «  YY, SEVAL) 

IMPLICIT  REAL*8  (A-H.U-Z) 

INTEGER  N 

DIMENSION  YIN).YY(N) 

C 

INTEGER  I , J, K 
DATA  1/1/ 

IFU.GE.N)  1*1 
IFCU.LT.Y(I))  GO  TO  10 
IF(U.LE.Y(I+1))  GO  TO  30 
10  1*1 

J*N*  1 

20  K*(I*J)/2 

IF(U.CT.Y(K))  J*K 
IFCU.GE.YCK))  I*K 
IF(J.i'T,I*l)  GO  TO  20 
C 

30  DY*U-Y(l) 

DLLYY*YY(I*1)-YY(I) 

DELY*Y(I*l)-Y(I) 

SEVAL»YY ( I ) *DY*DELYY/DELY 

RETURN 

END 

C  END  OF  SEVAL. 

SUBROUTINE  PD(W1,W2» TEMP, SMI 2 »0M11,PL1,PD2) 

C  PD  CAC'JLATES  PD (PRESSUF.E*01FFUsIUN  CONSTANT)  FOR  BINARY 

C  GAS  MIXTURE 

IMPLICIT  REAL*«(A-H,0-Z) 

PD 1*0. 002628* (TEMP**3«*(W1*W2) /(2,*hl*«2))**0.5/(SM12**2.0*OMll) 
PD2»PD1 
RETURN 
END 

C  END  OF  PD 
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